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» ■ We still do not know what causes aspherical planetary nebula morphologies. A plausible 

^ | hypothesis is that they are due to the presence of a close stellar or substellar companion. 

So far, only ^40 binary central stars of planetary nebula have been detected, almost all of 
them with such short periods that their binarity is revealed by photometric variability. Here 
we have endeavoured to discover binary central stars at any separation, thus determining the 
unbiased binary fraction of central stars of planetary nebula. This number, when compared 
to the binary fraction of the presumed parent population can give a first handle on the origin 
of planetary nebulae. By detecting the central stars in the I band we have searched for cool 

■ companions. We have found that 30% of our sample have an / band excess detected between 

one and a few a, possibly denoting companions brighter than M3-4V and with separations 
smaller than ^1000 AU. By accounting for the undetectable companions, we determine a 
de-biased binary fraction of 67-78% for all companions at all separations. We compare this 
number to a main sequence binary fraction of (50±4)% determined for spectral types F6V- 
G2V, appropriate if the progenitors of today's PN central star population is indeed the F6V- 
G2V stars. The error on our estimate cannot be constrained tightly, but we determine it to be 

v^a , between 10 and 30%. We conclude that the central star binary fraction may be larger than 

■ expected from the putative parent population. However, this result is based on a sample of 
' 27 bona fide central stars and should be considered preliminary. The success of the I band 

method rests critically on high precision photometry and a reasonably large sample. From a 
similar analysis, using the more sensitive J band of a subset of 1 1 central stars, the binary 
fraction is 54% for companions brighter than ^M5-6V and with separations smaller than 
about 900 AU. De-biassing this number in the same way as was done for the / band we obtain 
a binary fraction of 100-107%. The two numbers should be the same and the discrepancy is 
likely due to small number statistics. Finally, we note how the previously-derived short period 
PN binary fraction of 15-20% is far larger than expected based on the main sequence binary 
fraction and period distribution. 

As a byproduct of our analysis we present an accurately vetted compilation of observed 
main sequence star magnitudes, colours and masses, which can serve as a reference for future 
studies. We also present synthetic colours of hot stars as a function of temperature (20-170kK) 
and gravity (log g = 6 — 8) for Solar and PG1 159 compositions. 

Key words: planetary nebulae: general - binaries: general - stars: evolution - white dwarfs 
- techniques: photometric. 
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1 INTRODUCTION 

A single star may be incapable of generating non-spherical plane- 
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PN shapes (e.g. lGarcfa- SeguraetalJl999ll2005h have traditionally 
assumed the constancy of magnetic fields over the high mass-loss 
period (known as the superwind phase) that characterises the end 
of the asymptotic giant branch (AGB) evolution. However, even a 
weak magnetic field during the end of the AGB can act to slow 
down the differential rotation that generates the field in the first 
place: today we have no viable theory to sustain a magnetic field 
during the superwind phase in a sing le AGB star dSokeri l2006; 
iNordhaus et al .1 120071) . An alternative theory, that a binary com- 
panion might be responsible for the shaping action, has become 
central in the study of PN and formed the core of the binary hy- 
pothesis which postulates that PN form more readily around bi- 
naries, where by binary we mean a star accompanied by another 
star, a brown dw arf, or even a planetary system. For a review, see 
iDe MarcdJ2009h . 

There are several ways in which binary companions as light 
as planets can alter the shape of the AGB s uperwind and the 
subsequent PN dMastrodemos & Morris! 1 19991 : lEdgar et al]|2008l : 
IPassv et akl 12013) . However, this does not p rove that all non- 
spherical PN (about 80% of the entire sample; IParker et ai]|2006l) 
derive from a binary interaction. In order to determine the impact of 
binarity on PN formation a first, fundamental step is to determine 
the binary fraction of central stars of PN. 

Binary detection methods for central stars of PN have centred 
on the light variability technique, where a central, unresolved bi- 
nary, undergoes eclipses, suffers ellipsoidal distortion or where the 
cool companion is irradiated by the luminous hot one. This tech- 
nique is responsib le for the de t ection of almost all th e known cen- 
tral star binaries teondll2000l : iMiszalski et alj|2009r) . The binary 
fraction determined in this way is ~ 15-20%. This technique is 
biased against binaries with periods smaller than about 2 weeks, 
against binaries with the orbital plane near the plane of the s ky, and 
against companions with small radii dDe Marco et al . 2008). 

To determine the binary fraction for binaries with any or- 
bital separation we need a technique that is free of separation bi- 
ases, such as the detection of red an d infrared excess from pho- 
tometry or spectrosco py. So far only IZuckerman et alj d 199 ll) and 
Frew & Parker! J2007I) have carried out such studies - but see also 



Bentlevl Jl989h and lHolberg & MagargaJ fc005h . lZuckerman et all 
(1991) detected definitive A"-band excess in 50% of 30 central stars 
but concluded that only in three cases this could be ascribed to a 
companion, while in the others the emission may be due to hot 
dust. Frew & Parker (2007) analysed 32 objects with 2-micron All 
Sky Survey (2MASS) or Deep Near Infrared Survey of the South- 
ern Sky (DENIS) near-IR photometry and deduced that >53% of 
PN have a cool companion (the completeness limit of that survey 
and the error limits were not quantified). 

Key to the success of such survey are (i) extremely accurate 
photometry with well quantified uncertainties in at least two blue 
colours (e.g., B and V, to determine the reddening) and in at least 
one red colour (/ or J, to detect the excess flux due to the compan- 
ion); (ii) a sufficiently large sample (>150 objects) covering the 
majority of a volume-limited sample and (iii) the use of a red band 
that is not contaminated by dust, practically leaving only the I or 
J bands as feasible. Here we present the first paper in a series that 
uses / band photometry (and J where possible) to detect a red or 
near-IR excess. The / band is not as sensitive as the J band, but, 
if photometric errors can be limited to 1%, this method can be a 
very practical way to observe a substantial number of objects in a 
relatively short time. 

Ideally, we would like to derive the binary fraction and period 
distribution in separate samples, such as the non-spherical and the 



spherical subsets of PN. The need for a larger sample as well as 
the lack of a clear definition of morphological classes lead us to 
simplify this test: we simply derive the binary fraction (by which, 
hereafter, we intend the stellar binary fraction) among a volume 
limited sample of central stars of PN and compare it with that of 
the presumed main sequence progenitor population. The expecta- 
tion from the current scenario, whereby PN derive from stars in the 
mass range ~l-8 Mq, whether they are single or in binaries, is that 
the PN binary fraction should be slightly smaller than that of in- 
termediate mass main sequence stars, if companions down to the 
brown dwarf limit and at all separations can be sampled. A small 
difference between the binary population on the main sequence and 
during the PN phase is justified by the expectation that very close 
binary main sequence stars suffer a strong interaction on the red gi- 
ant branch (RGB) and do not ascend the AGB; also, it is expected 
that some binary interactions on the RGB or AGB could result in 
mergers. 

The binary fraction of the pr ogenitor population h as recently 
been measured to be (50±4)% (Ra ghavan et afl l201Cl) . This in- 
cludes all companions down to the planetary regime out to all sep- 
arations with primary stars in the spectral range F6V to G2V. We 
therefore expect a PN binary fraction of approximately 10 points 
lower. Clearly, if we are to explain 80% of all PN (those with non- 
spherical shapes) with a binary interaction, either the binary frac- 
tion of PN is higher than for the main sequence, and close to 80%, 
or, if it is not, a fraction of PN should be explained by planetary in- 
teractions. While testing for the presence of planets around central 
stars is not generally within reach yet, we start here with determin- 
ing the stellar-companion, PN binary fraction. 

In §[2]we describe our sample. In §[3]we present our observa- 
tions and data reduction, while in §[4] we present our measurements 
of the photometric magnitudes and their uncertainties. In § [5] we 
give the details of the technique to detect / and J band excess flux 
and the predicted accuracies and biases that derive from it. In § [6] 
we report our results, including a discussion of objects that were 
detected to be variable in the course of our observations. A com- 
parison between our results and the prediction for the single and 
binary scenarios follows in §[7] In §[8] we report details of individ- 
ual objects. We summarise and conclude in §|9] 



2 THE SAMPLE 

The goal of this survey is to determine accurate B, V and / band 
photometry of the ~200 central stars closest to the Sun. This sam- 
ple has been recently compiled using the best available data from 
the literature supplemented by an improved Ha surface brightness- 
radius relation which allows one to obtain distances accurate to 
~20% in most cases (Frew 2008). Obtaining the J band magni- 
tude of the central stars is also desirable but much less practical for 
such a large sample. 

The sample presented here consists of 30 central stars of PN 
which were selected solely based on their low PN surface bright- 
ness (radius of the PN is larger than ~25 arcsec in most cases) as 
well as on the faint V magnitudes of their central stars. The first 
criterion allows us to reduce the measurement error; the second in- 
sures that fainter companions can be detected because, although 
faint V magnitudes can imply intrinsically bright, distant objects, 
the large PN radius tends to select for closer objects whose faint V 
brightness indicates that the stars are intrinsically faint. 

In Tables Q] and 2 we report the best available literature data, 
as well as values determined here. We used trigonometric dis- 
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Name 


Sp. Type 


PN 

morph. 1 


D 

(kpc) 


My 

(mag) 


E(B - V) 


T eB (method 2 ) 
(K) 


logs' 5 


Reference & Comment 


A7 


DAO 


R? 


0.53 


6.8 


0.00 


99±18 (m) 


7.03±0.43 


Napiwotzki 1999; wide binary 


A 16 




R? 


2.50 


6.3 


0.11 


95±8 (z He II) 


5.1 or 7.5 




A 20 


0(H) 


R? 


2.35 


4.3 


0.09 


119±22(m) 


6.13±0.13 


Rauchetal. 1999 


A 28 


O(H) 


R? 


1.10 


6.2 


0.00 


70±30 (z He II) 


5.1 or 7.5 




A31 


DAO 


ISM 


0.62 


6.5 


0.04 


95±8 (z He II,m) 


6.63±0.3 


Napiwotzki 1999 


A 57 


O(H)" 


E/B 


3.0: 


4.0: 


0.38 6 


>60 (z H I) 


5.0 


"Miszalski etal. 2011 
6 Schlaflv & Finkbeiner 201 1 


A71 




ISM 


1.22 


7.7 


(0.8) 


145±10(zHen) 


6.2 or 7.0 




A 72 


PG1159 C 


E 


1.75 


4.8 


0.09 


> 100 (He II z) 


5.3 or 7.4 


c Douchin et al., in prep. 


A 79 


FO V 


T/B 


3.30 


2.5 


0.62 


165±25 (Tc) 


7.3 


F0V central star 


A 84 




E 


1.50 


7.3 


0.16 


101±10 (zHen.Tc) 


5.4 or 7.4 




DeHt5 


DAO 


- 


0.35 


7.5 


0.16 


70±10 (m) 


6.65±0.19 


Napiwotzki 1999; PN mimic 


EGB 1 


DA 


ISM? 


0.65 


6.7 


0.44 6 


147±25 (m) 


7.34±0.31 


Napiwotzki 1999: PN mimic? 


EGB 6 


DAO 


ISM 


0.60 


7.0 


0.03 6 


110±10 (m) 


7.36 


Gianninas et al. 2010 


HaWe5 


DA 


- 


0.42 


8.9 


0.20 


38±2 (m) 


7.58±0.20 


Napiwotzki 1999; PN mimic 


HDW3 


DAO 


ISM 


0.80 


7.0 


0.23 


125 ±28 (m) 


6.75±0.32 


Napiwotzki 1999 


HDW4 


DA 


- 


0.21 


9.6 


0.08 


47±2 (m) 


7.93±0.16 


Napiwotzki 1999; PN mimic 


IsWe 1 


PG1159 


ISM 


0.63 


7.0 




100±20 (m) 


7.0 


Werner 1995 


Is We 2 


DA 


E 


0.95 


7.3 




148±20 (z H I) 


5.4 or 7.4 




JnErl 


PG1159 


E/B 


1.15 


6.8 


0.06 


130±15 (m) 


7.0 


Rauch& Werner 09951 


K 1-13 




E/B 


2.45 


6.4 


0.03 


80±30 (z He II) 


5.4 or 7.4 


K2-2 


hgO(H) 


ISM? 


0.80 


4.7 


0.03 


69±15 (m) 


6.09±0.24 


Napiwotzki 1999; PN mimic? 


NGC 3587 


hgO(H) 


E 


0.78 


6.3 


0.00 


95±8 (m) 


6.94±0.31 


Napiwotzki 1999 


NGC 6720 


DAO 


E/B 


0.74 


6.3 


0.08 


125±25 (m) 


6.88±0.26 


Napiwotzki 1999 


NGC 6853 


DAO 


E/B 


0.41 


5.9 


0.07 


130±10 (m) 


6.72±0.23 


Napiwotzki 1999 


PuWe 1 


DAO 


R/B? 


0.37 


7.4 


0.09 


99±10(m) 


7.09±0.24 


see §HJ8] 


Sh 2-78 


PG1159 


E/B 


0.75 


7.3 


0.48 


120±15(m) 


7.4 


Dreizler 1999 


Sh 2-176 


DA 


ISM 


1.0: 


7.9 




150±25 (z He II) 


6.8 


Gianninas et al. 2010 


Sh 2-188 


DAO 


ISM 


0.83 


6.9 


0.28 


130±30 (m) 


6.82±0.6 


Napiwotzki 1999 


Ton 320 


DAO 


ISM 


0.55 


7.0 


0.00 


78±15 (m) 


7.76 


see §HTTJ 


WeDe 1 


DA 


ISM 


0.79 


7.5 




141±19 (m) 


7.53±0.32 


Napiwotzki 1999 



1 Legend: R: round; E: elliptical; B: bipolar; T: torus; ISM: features dominated by interaction with the inters tellar medium. 

2 Legend: m: spectral model; z He II: Zanstra method on helium lines; Z H I: Zanstra method on hydrogen lines; Tc: cross-over 
(Ambartsumyan) temperature. 

3 Values in Roman font are from stellar atmosphere models, while those in italics are derived from stellar evolutionary tracks, where 
we have chosen the higher of the two values for our calculations. 



Table 1. Observational data discussed in §§[2]and[8] 



Name J H K Data 

(mag) (mag) (mag) source 



A7 


16.10±0.08 


16.16±0.18 


> 15.08 


2MASS 


A31 


15.95±0.01 


15.80±0.0I 


15.67±0.01 


UKIDSS 


A 72 


16.67±0.13 






2MASS 


A 79 


15.05±0.04 


14.64±0.06 


14.44±0.08 


2MASS 


DeHt5 


15.57±0.07 


15.96±0.20 


15.58±0.22 


2MASS 


EGB 1 


16.64±0.16 


> 17.47 


> 15.62 


2MASS 


EGB 6 


16.52±0.10 


15.95±0.16 


16.10±0.26 


2MASS 




16.43±0.20 


16.08±0.09 


15.63±0.04 


FL93 


HaWe5 






17.78±0.18 


UKIDSS 


HDW3 




17.62±0.08 


17.49±0.12 


UKIDSS 


K2-2 


14.94±0.05 


14.99±0.06 


15.09±0.14 


2MASS 


NGC 3587 


16.71±0.13 






2MASS 


NGC 6720 


16.40±0.20 






2MASS 


NGC 6853 


14.75±0.05 


14.70 


14.61 


2MASS 


Sh 2-78 






17.89±0.15 


UKIDSS 


Ton 320 


16.59±0.16 


>15.96 


15.73±0.20 


2MASS 



Table 2. Near IR magnitudes from the 2-micron All Sky Survey 
(2MASS) and the UKIR T Infrared Deep Sky Survey (UKIDSS) and 
iFulbright & Liebenl Jl993l FL93) 



tances if available (e.g. Harris et al. 2007; Benedict et al. 2009), 
otherwise distances were taken from Frew (2008) and Frew et al. 
(2012). From the distances and de-reddened V magnitudes we de- 



rived the absolute magnitude, Mv ■ The effective temperatures are 
those derived via atmospheric spectrum fitting by several authors. 
If a model was not available we adopted the temperature provided 
by Zanstra analyses of either the hydrogen or helium lines. To do 
this we used the new V photometry reported here in combination 
with new integrated Ha fluxes taken from Frew (2008) and Frew et 
al. (2012). Since these only provide a lower limit to the temperature 
for optically-thin PN, we have used additional information, where 
appropriate, to determine the most suitable temperature value. The 
reddening values reported in Table Q] are derived from data in the 
literature other than the stellar B — V colour, i.e. from the nebular 
B aimer decrement or using the interstellar hydrogen column den- 
sity. In §[5] we will derive the reddening using a comparison of the 
observed B — V colour and that predicted for a single star of the ap- 
propriate temperature. At that time the nebular and stellar-derived 
values will be compared. 

Out of the 30 objects studied in the I band, three were 
later discovered to be PN mimics, nebulae which we exclude to 
have been ejected by the central star when it was on the AGB 
dFrew & Parker! l2010h . These are HaWe5 and HDW4, both of 
which have a very high surface gravity indicating a long cooling 
age. Such stars departed from the AGB well over the maximum life- 
time of a PN of ~ 100,000 years. Their nebulae must be Stromgren 
spheres ionised by the hot central stars. In addition we also consider 
DeHt 5 a PN mimic due to a range of discrepancies which we de- 
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Name 


RA 


Dec 


Exp Times 
B V R I 
(sec) 


UT of the observation 


Night 


Comment 




ns m m o 

Uj U j U / jz 


ic 1A 1~) S 


fin fin fin 

OU, OU, DU, — 

45,45.45,90 

AC\ Aft AC\ sn 
4U, 4U, 4U, OU 


1 l.ZZ JU/1U/ZUU/ 

08:19 03/11/2007 
ns-is ns/i i nnm 

US.JJ U j/ J J/zUU/ 


5 


No / band 


A 16 


06 43 55.46 


+61 47 24.7 


120, 120, 120, 240 
150, 150, 150, 300 

1 lf\ 1 1C\ 1 on 1A A 

JZU, JZU, IZU, Z4U 


09:56 02/1 1/2007 
07:50 03/1 1/2007 

i n- ^ o A/i / 1 i /inm 
JU.jy U4/I J/zUU/ 


4 
5 
6 




A 20 


07 22 57.74 


+01 45 32.8 


120, 120, 120, 240 
120, 120, 120, 120, 240 
i nn i nn i nn 9nn 

JUU, 1UU, 1UU, ZUU 


10:49 0271 1/2007 

09:02 03/1 1/2007 

no- 19 ns/ 1 i /9nm 
uy. JZ Uj/ j j/zuu / 


4 
5 




A 28 


08 41 35.57 


+58 13 48.4 


150, 150, 150, 300 
120, 120, 120, 240 

i in i ia ion iAt\ 
IZU, JzU, JzU, Z4U 


10:15 03/11/2007 
11:59 04/11/2007 

i i - i q / 1 i /oAm 
J J.Jft Uj/J J/zUU/ 


5 
6 




A31 


08 54 13.16 


+08 53 53.0 


45,45,45,90 
40, 40, 40, 80»2 


11:24 03/11/2007 
10:44 05/11/2007 


5 
7 




A 57 


10 IT IIS 7) 

i y i / u j . i j 


+25 37 33 4 


i nn i nn i nn i sn 

1UU, JUU, JUU, 1JU 

1 9n 1 9n 1 9n iac\ 

IZU, IzU, JZU, ZH-U 


n^-m n9/i i/9nm 
jjj.uz uz/ j j/zuu/ 

n^-A9 n^/i i /9nm 

UJ.+Z Uj/ j J/ZUU / 






A71 


20 32 23.22 


+47 20 50.4 


150, 150, 150, 300 
9nn 9nn 9nn 9nn 

ZUU, ZUU, ZUU, ZUU 

250, 250, 250, 500 


02:52 03/11/2007 
nA-99 nzj/i i nam 

U+.ZZ U+/ J J/ZUU / 

02:40 05/1 1/2007 


5 
7 




A 72 


?n sn n9 nfi 

ZU JU uz.uo 


4.1 1 11 9Q fi 
tij zy.o 


sn sn sn sn 

OU,OU,OU,OU 

60, 60, 60, 120 
fin fin fin i 9n 

OU, OU, OU, JzU 


n^-9d ni /m /9nm 

UJ.ZH U 1/ 1 17 ZUU / 

04:59 02/11/2007 
c\a-ii n*>/i i/9nm 

JJ+.JJ U j/ J J/zUU/ 


4 




A 79 


22 26 17.27 


+54 49 38.2 


100, 100, 100, 200 
100, 100, 100, 200 
9nn 9nn ?nn isn 

ZUU, ZUU, ZUU, JJU 


04:14 01/11/2007 
05:35 02/1 1/2007 
nfi- i q nzi/i i nnm 

UO. 1 y UM7 1 1 / ZUU / 


3 
4 




A 84 


11 Al AA (Y) 


i C 1 11 S£ Q 

+j i Zj jo.y 


sn sn sn i nn 

oU, oU, oU, JUU 

80, 80, 80, 160 
120, 120, 120, 240 


nA-^ i n i / ii nnm 

JJ+. J J U il J J/zUU / 

06:02 02/1 1/2007 
05:24 05/1 1/2007 


4 
7 




vrn f 


ZZ IV JJ. / J 


+ /u jo uj. i 


cn cn sn sn 

JU, jU, jU, oU 


n^-^9 n9/i i nnm 
UJ. JZ uz/ j j/zuu / 








n i mm so 
ui u / u / .jy 


+ / J J J Zj.Z 


sn sn sn i fin 

ml, oU, OU, JOU 

sn sn sn i fin 

ml, oU, oU, JOU 


nfi-^zj n9/ 1 i /9nm 

UO.J4 UZ/ J J/ZUU/ 

UJ.JZ u j/ j j/zuu/ 






EGB 6 


no S9 SS QQ 
uy jz jo.yy 




fin fin fin i 9n 

OU, OU, OU, IZU 

fin fin fin i 9n 

OU, OU, OU, JZU 


1 9- ii n9/i i /9nm 

1 Z. 1 J V)z./ I l/ZUU / 

i i -no n^/ 1 i nnm 
J i .uy UJ/ 1 J/ZUU / 






HaWe 5 


UJ 4j Z0.04Z 


+3 / 48 J J . / 


1 AA 1 AA 1 AA OflA 

JUU, JUU, JUU, ZUU 

100, 100 100. 200 


AQ-/1A AT / 1 1 11 Am 

U8.4U UZ/J J/ZUU/ 
05:34 03/1 1/2007 


5 




HDW 3 


97 IS Ad 


j-zlS 9zL 90 S 

T^J Z4- ZU.J 


1 9n i 9n 1 9n 9zin 

1 ZU, 1 ZU, 1ZU ? ZM-U 

160, 160, 160, 320 
100, 100, 100, 200 


n7-^.n in/m/9nm 

\) 1 .H-VI JU/H^/ZUU/ 

08:01 02/11/2007 
06:47 05/1 1/2007 


j 

4 

7 






AS 17 Sfi 11 
Uj j / jo.ZJ 


i« 19 i fi n 

TJJ JZ 10. u 


1 9n 1 9n 1 9n 9zin 

JZU, JZU, JZU, Z4U 


no-^n n9/ 1 i nnm 
uy.ju uz/ j j/zuu/ 






IsWe 1 


03 49 05.89 


+50 00 14.8 


120, 120, 120, 240 
100, 100, 100, 200 

1 AA 1 AA 1 AA 1AA 

JUU, JUU, JUU, ZUU 


05:29 03/1 1/2007 
06:58 04/1 1/2007 

A"? - ii a<;/ i i ii Am 
U /. JZ UJ/ J J/ZUU / 


5 
6 




IsWe2 


22 13 22.53 


+65 53 55.5 


80, 80, 80, 160 


03:48 01/1 1/2007 


3 










80, 80, 80, 160 


04:13 02/11/2007 


4 


/ band done later at 4:53 








sn sn sn i fin 

5U, oU, oU, JOU 


m-9^ nin i nnm 

UJ.ZJ U J/ J J/ZUU/ 












150, 150, 150, 250 


05:47 04/1 1/2007 


6 










i nn i nn i nn 9nn 

JUU, JUU, JUU, ZUU 


nA-^on^/i mnm 
U4.jy uj/ j j/zuu/ 






JnEr 1 


07 57 51.63 


+53 25 17.0 


120, 120, 120, 120 
100, 100, 100, 200 


11:20 02/11/2007 
09:43 05/11/2007 


4 
7 






AO f\& AA £A 

US 06 46. 3U 


m £1 7/1 Q 
-UZ 32 34.8 


i <n i cn i ca iz.c\ 
IjU, JjU, 1jU, ZjU 

120, 120, 120, 200 
i sn i sn i sn 9sn 

JjU, JjU, JjU, ZjU 


i i - «; /i a i / 1 i/i Am 
1 J.J4 UJ/J J/ZUU/ 

11:43 02/11/2007 
1 1 . «n a^5/ i i nnm 

J J . JU UJ/ J J/ZUU / 


3 
4 




K2-2 


06 52 23.17 


+09 57 55.7 


30, 30, 30, 60 
ic\ in ic\ fin 

JU, JU, JU, OU 

30, 30, 30, 60 


10:29 02/1 1/2007 
nfi-Ai n^/i i nnn7 

Uo.M-l Uj/ L l/ZUU/ 

08:57 05/11/2007 


4 
7 




NGC 3587 


11 14 47.73 


+55 01 08.5 


60, 60, 60, 120 
fin fin fin i 9n 

OU, OU, DU, JzU 

60, 60, 60, 120 


12:35 02/11/2007 
1 9-97 c\ai i i nnm 

JZ.Z / U4/ J J/ZUU / 

12:41 05/11/2007 


4 
6 




NGC 6720 


18 53 35.08 


+33 01 45.0 


50, 50, 50, 50 
in if\ in fin 

JU, JU, JU, OU 


01:59 02/11/2007 
m-9i n*>/i innn? 

UJ.ZJ UJ/ J J/ZUU/ 


4 




NGC 6853 


19 59 36.38 


+22 43 15.7 


15, 15, 15, 15 
in if\ if\ fin 

JU, JU, JU,OU 

20,20,20,40 


03:32 02/1 1/2007 
n9-^zj mi i i nnn7 

UZ. J4 UJ/ J J/ZUU / 

04:13 05/11/2007 


4 
7 




PuWe 1 


06 19 34.33 


+55 36 42.3 


50, 50, 50, 100 
50, 50, 50, 100 


07:25 03/1 1/2007 
07:40 05/11/2007 


5 

7 




Sh 2-78 


19 03 10.09 


+14 06 58.9 


100, 100, 100, 150 
150, 120, 120, 400 


02:32 02/11/2007 
01:58 03/11/2007 


4 
5 




Sh 2-176 


00 31 54 


+57 22.6 


80, 80, 80, 160 
80, 80, 80, 160 
150, 150, 150, 300 


05:17 01/11/2007 
06:27 02/1 1/2007 
05:50 05/1 1/2007 


3 
4 
7 




Sh 2-188 


01 30 33.11 


+58 24 50.7 


120, 120, 120, 240 
120, 120, 120, 240 
120, 120, 120, 240 


07:16 02/11/2007 
04:17 03/11/2007 
06:21 05/11/2007 


4 
5 
7 




Ton 320 


08 27 05.53 


+31 30 08.6 


60, 60, 60, 120 
60, 60, 60, 180 


09:42 03/11/2007 
10:15 05/11/2007 


5 

7 




WeDe 1 


05 59 24.87 


+10 41 40.4 


120, 120, 120, 240 
120, 120, 120, 240 
120, 120, 120, 240 


06:53 03/1 1/2007 
11:31 04/11/2007 
08:05 05/1 1/2007 


5 
6 
7 





Table 3. Log of observations taken in photometric conditions 
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scribe fully in § 18.81 Although we suspect that two other objects are 
also PN mimics (EGB 1 and and K 2-2), for the moment we keep 
them in the central star sample. Our bona fide sample comprises 27 
objects studied in the I band of which 1 1 have additional J-ban d 
photometry. A 79 has a cool spectral type ( Rodri guez et"ai1l200ll) . 
indicating that it is a binary, but this was not realised at the time of 
the observations. 

In Column 3 of Table[T]we have listed the morphologies of the 
PN. Many of them are too old and dispersed for a morphology to be 
determined accurately, and the PN features are dominated by inter- 
action with the ISM. Those PN we list as "R?" are approximately 
round, but may contain internal structure which, in some cases, is 
distinctly bipolar. An example is A 7, which, although diffuse, con- 
tains a hint of bipolarity in the form of brighter structure on either 
side of the st ar. None of the PN we mark as "R?" is like the round 
PN Abell 39 Jjacobv et al]200lh . A 20 is a faint round PN with an 
evacuated cavity around the central star. Its appearance is that of 
a pole-on doughnut. A 28 on the other hand, has a sharp, broken 
rim and may be more elliptical than round. We will not mention 
morphology further. It is not good form, when trying to determine 
if binarity causes certain morphologies, to use the morphology as 
an argument for binarity. 



3 OBSERVATIONS AND DATA REDUCTION 

The observations were acquired during 8 nights between the 30th 
of October and the 6th of November, 2007 at the 2.1-m tele- 
scope at Kitt Peak National Observatory. However, the data from 
nights 2 and 8 were not photometric. The weather conditions dur- 
ing the other nights were mostly photometric. Data taken in non- 
photometric conditions were used to monitor stars for relative vari- 
ability (see § 16. lb . A log of the observations is presented in Table[3] 
The detector was a 2048 x 2048 pixel Tetronix CCD, which 
was binned 2 x 2 (to achieve a faster read-out and reduce the read 
noise level). The pixel size was 24 /im, and the field of view was 
10.2' x 10.2' (the platescale is 0.60 arcsec per (binned) pixel). The 
electronic gain of the camera was 3. 1 e~ / ADU (Analog-to-Digital 
Unit), which minimally samples the system noise of 6 e - , while 
providing a maximum signal of about 200,000 e~ before saturat- 
ing. 

The observations were made through the B,V,R and / 
Johnson-Cousins astronomical bandpasses. Ten bias frames per 
nights were obtained as well as ten dome flats per filter. Their me- 
dians were used to de-bias and flat fi eld the expos ures. Standard 
stars were selected from the list of lLandoltl Jl992h . so as to en- 
compass the colours and brightnesses of our targets as well as to 
sample a range of air masses. They are listed in Table 4 alongside 
with the nights when they were used. We applied a shutter correc- 
tion of —0.001 sec to the header exposure times, very low for even 
our shortest exposure times. 



4 THE DETERMINATION OF THE PHOTOMETRIC 
MAGNITUDES 

Instrumental magnitudes were m easured with the APPHOT pack- 
age in IRAl(3 hbdvll 198611 19931) . The radius of the aperture that 
samples the flux was estimated to be 8 binned pixels by measuring 



bright stars through increasingly larger apertures. The background 
was sampled between 8 and 13 binned pixels from the position of 
the star. Targets were measured using apertures of 3 binned pixels 
and then aperture-corrected. However, it was noticed that the final 
result was not very different from that obtained by measuring tar- 
gets through an 8-binned pixel aperture directly with no aperture 
correction. 

The instrumental magnitudes b,v,r and i, were then con- 
verted to observed magnitudes, B,V,R and I on the standard sys- 
tem, by making the following first order transformation: 

B = Ob + b + C B (B - V) - K B * Z B 
V = Ov + v + C v \b-V)-K v *Zv 
R = R + r + C R (V-R)-K R *Z R 
/ = Oi + i + d(R -I)-Ki*Zi 

where Ob,Ov,O r ,Oi are the instrumental offsets, 
C B ,Cv ,Cn,d are the colour terms, Kb,Kv,K r ,Ki are 
the extinction coefficients, and Z is the airmass (where the Z 
values are all the same and the subscripts are there for clarity). 
There can be a second order correction due to the atmosphere, 
and therefore zenith distance because the atmosphere acts as a 
broadband colour filter. This second order correction was negligi- 
ble. Although the airmass of our observations varied, it was never 
higher than 2. 1 and generally very close to unity. 

Our standard stars were used to calculate the needed coef- 
ficients. Each standard star was observed through each filter at 
several values of the airmass, Z, and its instrumental magnitudes 
b, v, r and i were determined. Then, for each filter, we solved the 
system of equations in QJ using a least squares method, where the 
unknowns are K, C and O. We finally plot Y\ = V — Ov — v — 
Gy (B — V) as a function of Z, and Y2 — I — Or — i + Ki * Z as a 
function of R— I. In Figure[T] we show the importance of obtaining 
observations of several standards at different airmass values. When 
there are few stars (less than ~25 stars), the fits tend to be notice- 
ably poorer. The values obtained for each coefficient as well as the 
total number of stars in each filter are summarised in Table 5. 



4.1 Uncertainties 

Several precautions are taken to minimise the uncertainty. First, 
the targets' signal-to-noise ratio is >100. Second, multiple obser- 
vations of each target are taken on different nights. This guards 
against possible mistakes in estimating the photometric conditions. 
It also allows us to assess variability affecting 15-20% of all cen- 
tral stars of PN and usually denoting binarity in its own right (see 
^ I6.lt . Finally, most of our targets have faint, large PN (>25 arc- 
sec), which reduces the error due to background subtraction. 

The uncertainties on the instrumental magnitude measure- 
ments (<Tb, a v , o r , and Oi) provided by IRAF include the photon 
statistics and the uncertainty on the sky measurement. While the 
shutter correction was included, we did not include a shutter cor- 
rection error because it was deemed negligible. By adding the un- 
certainties in quadrature we have, e.g., in the case of the error on 
the B-band magnitude: 



IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 



Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Coeff 


Ni^ht 1 


Night 3 


Night 4 


Night 5 


Ni°ht 6 


Night 7 


Ob 


— 0.8638 ± 0.0083 


11 U '1(1(1 (1 (1 1 

— U.82U9 ± U.U157 


(1 O AQH 1 (1 (1(1U(1 

—0.8436 ± 0.0080 


—0.8525 ± 0.0088 


A (1(11 1 _L n (11 CO 

—0.9012 ± 0.0153 


(1 OICC _L f\ 1^17(1 

—0.8466 ± 0.0079 


Ov 


—0.8054 ± 0.0073 


—0.7767 ± 0.0168 


—0.7877 ± 0.0100 


—0.7958 ± 0.0069 


—0.8155 ± 0.0167 


—0.7970 ± 0.0064 


Or 


-0.7940 ±0.0068 


-0.6169 ± 0.0244 


-0.7656 ±0.0127 


-0.7824 ± 0.0140 


-0.7775 ± 0.0193 


-0.7887 ±0.0127 


0, 


-1.774 ±0.008 


-1.790 ± 0.038 


-1.756 ±0.013 


-1.778 ±0.017 


-1.758 ±0.035 


-1.752 ±0.013 


Kb 


0.2286 ± 0.0059 


0.2550 ± 0.0107 


0.2370 ± 0.0054 


0.2355 ± 0.0056 


0.2205 ± 0.0106 


0.2522 ± 0.0055 


K v 


0.1313 ±0.0050 


0.1445 ± 0.0113 


0.1373 ±0.0070 


0.1325 ± 0.0043 


0.1441 ±0.0101 


0.1464 ±0.0042 


Kr 


0.09011 ± 0.00468 


0.2310 ± 0.0165 


0.1011 ±0.0091 


0.08112 ± 0.00886 


0.1072 ± 0.0178 


0.09376 ± 0.00863 


K, 


0.04785 ± 0.00537 


0.05281 ± 0.02654 


0.04264 ± 0.00956 


0.01315 ±0.01101 


0.04623 ± 0.02283 


0.04895 ± 0.00935 


Cb 


0.02832 ± 0.00320 


0.03006 ± 0.00546 


0.02772 ± 0.00308 


0.03248 ± 0.00298 


0.02849 ± 0.00600 


0.02747 ± 0.00273 


Cv 


-0.03070 ± 0.00264 


-0.02954 ± 0.00607 


-0.03141 ±0.00335 


-0.03480 ± 0.00232 


-0.03465 ± 0.00670 


-0.03137 ± 0.00254 


Cr 


-0.03581 ± 0.00472 


-0.05941 ± 0.01434 


-0.04464 ± 0.00754 


-0.06885 ± 0.00696 


-0.1066 ±0.0125 


-0.06955 ± 0.00804 


c, 


0.05139 ± 0.00586 


-0.2063 ±0.0161 


0.03163 ± 0.00944 


-0.06926 ±0.01115 


-0.1347 ±0.0258 


-0.02935 ± 0.00943 


# Stars 


27,28,24,23 


21,18,11,9 


54, 47, 35, 27 


46, 49, 33, 35 


24,28,22,25 


50,63,54,50 


(B,V,R.I) 















Table 4. Coefficients and errors for every photometric night. As explained in §[3] nights 2 and 8 were non-photometric 



-0.25 
-0.3 

-0.35 
-0.4 

-0.45 
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Night 4, filter B 



1.2 



1.4 



1.6 



-0.08 
-0.1 

-0.12 
>-~" -0.14 

-0.16 
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-0.2 
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Night 4, filter R 



1.2 



1.4 1.6 
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0.15 
0.1 
0.05 


-0.05 



0.5 

B-V 



0.02 


-0.02 
-0.04 
-0.06 



2 



0.2 0.4 0.6 
V-R 
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Figure 1: Example of fits to standard stars to estimate conversion coefficients. Note the small range on the y-axis of the R band fit 



2 

o"s = 

°o B +ol + {B - V) 2 crc B + C 2 B o 2 B _ v + Z 2 b o 2 Kb + al p 



where <j„ 



0.004 mag is the error on the aperture correction, 



calculated as the standard deviation of the aperture correction for 
the 5 reference stars used to calculate the aperture correction itself. 
By adopting the approximation o 2 B _ v — o 2 B + <jy « 2<jy, we 
obtain the following equations for the errors on the B and V bands: 



{°l B + °l + (B — V) 2 a 2 CB + Z 2 B a 2 KB + a 2 pco ) / 
(1 - 2C%) 



(2) 



cr 2 v = 
2 



(°"CV + o- v + {B - V) a Cv + Z v a Kv + a apco ) / 
(1 - 2C 2 ) 



which we solve simultaneously to derive the values of a B and ov . 
Similarly, we derive the uncertainties on R and I. 

When more than one independent observation was taken for a 
given target, the weighted mean was calculated: 



/' : 



(3) 



where N is the number of observations, Xi is a given measurement 
and pi is its associated probability (p; = (l/(7i)/Ej (l/<jj), where 
Gi is the error on that measurement and Oj is the error on the j 
measurement). The error on the mean was calculated in the follow- 
ing way: 



E£ lPi (*i-M) 2 



(4) 



The photometric magnitudes and their errors thus obtained are re- 
ported in Table[6] along with the number of measurements that were 
used in obtaining these values. In Appendix A we report the indi- 
vidual values and their individual errors. 
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Name 


# stars 


Night observed 


GD246 


IB 


1,3,4,5,6,7 


PG023 1+051 


IB + 6 


1,3,4,5,6,7 


PG09 18+029 


IB + 4 


5 


PG1047+003 


IB + 3 


5 


PG22 13-006 


IB + 3 


3,4,5,6,7 


RU 149 


2B + 6 


1,4,5,6,7 


92 245 1 


4 


1,3,4,5,6,7 


95 112 2 


IB + 4 


1,3,4,5,6,7 


98 650 3 


IB + 9 


1,3,4,5,6,7 



1 Tbis field also contains 92 248, 92 249 and 92 250. 
2 This field also contains 95 41-43 and 95 1 15. 
3 This field also contains 98 653, 98 670-671, 
98 675-676, 98 L5, 98 682, 98 685 and 98 1087. 



Table 5. Standard stars from lLandolJ \ 19921) and lLandolJ 120091) . In Column 
2 we indicate the number of blue stars (marked "B"), as well as of redder 
stars 



5 BINARY DETECTION TECHNIQUE BY RED AND IR 
EXCESS FLUX 

The ideal spectral location for this technique is the J band. This is 
the best compromise between companion brightness and the elimi- 
nation of contamination from hot dust. However, procuring J band 
data has proven more challenging than optical data because IR in- 
struments are often associated with larger telescopes which tend to 
allocate shorter observing runs and because of the need for pho- 
tometric conditions. In addition, for most of our sample we also 
need to obtain B- and V-band photometry because the values in 
the literature tend not to be sufficiently accurate, thus generating 
the need for parallel proposals to more than one telescope. As a 
result, we have found it more practical to use the / band to detect 
companion-generated flux excess as our work horse, accompanied 
by the J band wh en available (we used UKIRT, Infrared Deep Sky 
Survey (UKIDSS uLawrence et al.ll2007l) data, or the most reliable 
2MASS values - see Table 2). Below we discuss the technique and 
its biases. 

Most of our targets have reasonably well determined effec- 
tive temperatures either via stellar spectrum modelling or using the 
Zanstra technique (Table [T). Those that were modelled also have 
known gravities. For those stars with only a Zanstra temperature 
estimate, we determine the gravity using the most common cen- 
tral star mass, namely O.6I-M0, sel ecting the correspo nding stel- 
lar evolution ary track from fig. 2 of Napiwotzki ( 199 9j, who used 
trac ks from ISchonbernetl dl983h . iKoester & Schonberned Jl986h 
an d lB16cked ( fl995h ) and choosing the larger of the possible gravity 
values, appropriate for our sample of evolved PN. Once the temper- 
atures are obtained, single star colours (B — V,V — I and V — J) 
can be predicted. The predicted B — V colour is used together 
with the measured one to determine the reddening, E(B — V). 
If the reddening thus determined has a negative value, possible be- 
cause of random errors, we set the value to zero. As we will see in 
§ [6] all negative reddening values we derive in this way are very 
small and within the uncertainty, reassuring us that the uncertain- 
ties have been assessed reasonably. With the reddenings thus de- 
rived we obtain valu es of A\/E(B — V) using the reddening law of 
ICardelli et al.l ( fl989l) , where A represent the bandpass central wave- 
lengths. The central wavelengths for the filter bandpasses are ob- 
tained by convolving the bandpasses with a synthetic stellar atmo- 
sphere with T c ff =100 kK, log g = 7 and solar abundance. This de- 
creases the central wavelengths by approximately 10 A, compared 
to the values obtained for un-convolved bandpasses, but yields only 



a very small change in the final results. The bandpass central wave- 
lengths thus determined and the values of the wavelength-specific 
extinctions are presented in Table|7] 

The V — I, or V — J observed dereddened colours are then 
compared with the predicted ones so as to determine if a flux ex- 
cess exists. An I (J) band excess is detected every time the ob- 
served (dereddened) and predicted colours are different by more 
than the combined uncertainties. While a high confidence result 
can already lend substantial weight towards a binary interpretation, 
a lower sigma result needs to be confirmed by additional photome- 
try or spectroscopy. 

To predict the B - V, V - I, V - J, R - I and J - H 
colours of single post-AGB stars, we use theoretical stellar atmo- 
sphere models calcula ted with the simulation code TMAW, the web 
interface to TMAP iWerner & Dreizlej|l999l ; IWerner et al]|2003l : 
iRauch & Deetienll2003h . or the German Astrophysical Virtual Ob- 
servatory grid calculations TheoSS^fl The colours and further de- 
tails of the calculations are provided in Appendix |B1 

The uncertainties on the measurements (see § 14.1b are com- 
bined with the uncertainty on the predicted colours to derive an 
uncertainty on the colour excess. The uncertainty on the theoretical 
colours reflect solely the uncertainty on the temperature. We as- 
sume that additional sources of uncertainty are far smaller and do 
not play a role. The uncertainty in colours derived from the uncer- 
tainty on the temperature is quite small, in particular for stars hotter 
than ~50 kK, which is the case for most of our sample. Systematic 
uncertainties on the t heoretical, single star colours are estimated 
to be below 1% (e.g- lRauch et al . 2007). In fact, the observational 
uncertainties dominate the error budget. 

To further exemplify the technique, we present in Fig. [2] 
the predicted I- and J-band excess as a function of primary 
star absolute V magnitude (My) and main sequence companion 
mass/spectral type. This plot gives a good idea of the ability of the 
I- and J-band methods to detect companions, although it must be 
pointed out that this figure is only valid for those stars that have 
already entered the cooling track. To generate Fig. [2] we created a 
grid of Mv absolute stellar brightness for the hot central stars using 
the effective temperatures to determine both bolometric luminosi- 
ties and bolometric corrections. The former was determined using 
th e average of the coo li ng track temperature - lumin osity relations 
of ISchonberned dl993h . IVassiliadis & Wood H994I) and iBlocked 
(1995) (for a 1.5 M© main sequence star, but the scatter is very 
small). The latter was determined us ing the average of the rela- 
tion determined bv lVacca et alj and a blackbody curve. To 
each hot central star we added the flux of a range of main sequence 
companions using magnitudes and colours reported in AppendixICl 
We thus created a grid of binaries. Next we calculated the total 
(binary) B — V colour for each primary-secondary combination 
and compared it to that of the primary alone. A difference is only 
present when the companion contributes flux in the V or even B 
and V bands, which is the case for companions brighter than K0- 
5 V. Any difference is interpreted as reddening and used to deredden 
all bands. We then subtract the primary I band flux from the binary 
(dereddened) one and any difference is labelled as excess (this is 
the value reported on each contour line). 

As can be seen in Fig. [2] a measured I or J band excess cor- 
responds to two distinct companion spectral types. This is due to 
the contribution of the brightest companions to the V and even B 
bands, resulting in too high a predicted reddening, too blue a de- 
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Name 


B 


V 




R 




I 




A7 


15. 190 ± 0.003 ( 3 ) 


15.495 ± 0.004 


3 ) 


15.632 ± 0.013 


3) 


15.818 ±0.011 


2) 


A 16 


18.517 ± 0.015(3) 


18.714 ±0.018 


:3) 


18.733 ±0.020 


3) 


18.686 ±0.016 


'3 ) 


A 20 


16.216 ± 0.002 ( 3 ) 


16.466 ± 0.006 


3 ) 


16.559 ± 0.007 


3) 


16.701 ±0.016 


;3> 


A 28 


16.280 ± 0.008 ( 3 ) 


16.557 ± 0.009 


3) 


16.691 ± 0.008 


3) 


16.877 ±0.014 


;3 ) 


A31 


15.201 ± 0.007 ( 2 ) 


15.544 ± 0.001 


2) 


15.693 ± 0.009 


2) 


15.831 ±0.016 


'2) 


A 57 


17.903 ± 0.001 ( 2 ) 


17.734 ±0.011 


,2) 


17.451 ±0.010 


2) 


17.210 ±0.002 


;2> 


A71 


19.384 ±0.015(3) 


19.335 ± 0.006 


3) 


19.253 ± 0.010 


3) 


19.185 ±0.056 


;3 ) 


A 72 


15.761 ±0.021 (3) 


16.070 ± 0.028 


3) 


16.237 ± 0.091 


3) 


16.381 ±0.043 


;3) 


A 79 


17.825 ±0.013 ( 3 ) 


16.965 ± 0.005 


,3) 


16.397 ± 0.026 


3) 


15.743 ± 0.085 


'3) 


A 84 


18.366 ±0.012(3) 


18.584 ±0.013 


3) 


18.613 ±0.020 


3) 


18.671 ±0.020 


'3 ) 


DeHt5 


15.268 ±0.018 ( 1 ) 


15.495 ± 0.019 


1) 


15.568 ±0.018 


1 ) 


15.631 ±0.018 


; i ) 


EGB 1 


16.308 ± 0.007 ( 2 ) 


16.439 ± 0.005 


2) 


16.452 ± 0.007 


2) 


16.482 ± 0.022 


;2> 


EGB 6 


15.692 ± 0.002 ( 2 ) 


15.999 ± 0.002 


2) 


16.137 ±0.008 


2) 


16.300 ± 0.009 


'2) 


HaWe5 


17.321 ±0.014(2) 


17.439 ± 0.005 


,2) 


17.471 ± 0.004 


2) 


17.528 ±0.011 


;2> 


HDW 3 


17.084 ± 0.003 ( 3 ) 


17.187 ±0.004 


3 ) 


17.218 ±0.039 


3) 


17.234 ±0.019 


;3 ) 


HDW4 


16.310 ±0.011 ( 1 ) 


16.540 ±0.013 


1) 


16.638 ±0.017 


1 ) 


16.739 ±0.017 


; i ) 


IsWe 1 


16.374 ±0.017(3) 


16.523 ± 0.007 


3) 


16.576 ± 0.016 


3) 


16.644 ±0.013 


'3 ) 


Is We 2 


18.142 ±0.026(5 ) 


18.160 ±0.033 


5) 


18.118 ±0.026 


5) 


18.098 ± 0.022 


;5 > 


JnEr 1 


16.775 ± 0.005 ( 2 ) 


17.128 ±0.013 


2) 


17.288 ± 0.001 


2) 


17.501 ± 0.023 


;2) 


K 1-13 


18.051 ±0.013 ( 3 ) 


18.425 ±0.006 


3) 


18.592 ±0.019 


3) 


18.846 ±0.044 


;3) 


K2-2 


13.977 ± 0.007 ( 3 ) 


14.263 ± 0.010 


:3> 


14.390 ± 0.008 


3) 


14.553 ±0.013 


'3 ) 


NGC 3587 


15.414 ±0.001 (3) 


15.777 ± 0.009 


13) 


15.960 ± 0.006 


3) 


16.194 ±0.029 


'3 ) 


NGC 6720 


15.405 ±0.016(2) 


15.769 ± 0.023 


2) 


15.901 ± 0.003 


2) 


16.062 ± 0.012 


;2) 


NGC 6853 


13.749 ± 0.026 ( 3 ) 


14.089 ± 0.010 


3 ) 


14.247 ± 0.006 


3) 


14.405 ± 0.010 


;3 > 


PuWe 1 


15.291 ± 0.008 ( 2 ) 


15.545 ± 0.006 


2) 


15.662 ±0.011 


2) 


15.792 ± 0.008 


;2> 


Sh 2-78 


17.633 ±0.012 (2) 


17.660 ± 0.005 


2) 


17.608 ± 0.025 


2) 


17.543 ± 0.027 


'2) 


Sh 2-176 


18.489 ±0.086( 3 ) 


18.559 ±0.019 


3) 


18.570 ±0.016 


3) 


18.545 ±0.029 


;3) 


Sh 2-176 1 


18.442 ±0.018 (2) 


18.551 ±0.016 


2) 


18.562 ±0.009 


2) 


18.531 ±0.010 


;2) 


Sh 2-188 


17.424 ± 0.013(3) 


17.447 ± 0.004 


3) 


17.398 ± 0.009 


3) 


17.376 ±0.011 


'3) 


Ton 320 


15.379 ± 0.007 ( 2 ) 


15.725 ± 0.006 


2) 


15.890 ± 0.010 


2) 


16.105 ±0.018 


'2) 


WeDe 1 


16.958 ± 0.007 ( 3 ) 


17.226 ± 0.004 


3 ) 


17.338 ±0.010 


3) 


17.489 ± 0.016 


3) 



These measurements are obtained by excluding the observations taken in night 3. 



Table 6. The photometric magnitudes of our targets rounded to 3 decimal places. Formal errors lower than 1% were set at 1%. In brackets are the number of 
independent exposures (taken on different nights) used to calculate the final photometric magnitude and the uncertainty. In Table lAll we present the individual 
magnitudes whose weighted averages are displayed here 



Band 


Ao 


A X /E(B - V) 


U 


3597 A 


4.86 


B 


4386 A 


4.12 


V 


5491 A 


3.10 


H 


6500 A 


2.10 


1 


7884 A 


1.90 


J 


1.237 fim 


0.889 


H 


1.645 fim 


0.562 


K 


2.212 fim 


0.349 


J2MASS 


1.241 fim 


0.885 


-H2MASS 


1.651 /mi 


0.349 


-ft"2MASS 


2.165 fim 


0.361 



Table 7. Bandpass central wavelengths after convolution with a 100 kK, 
logg=7, solar abundance synthet ic stel lar atmosphere and resulting extinc- 
tions according to Car delli et alj jl989t) 

reddened spectral energy distribution (SED) and therefore smaller 
I or J-band excess. This effect is less pronounced for the J band, 
because reddening effects are smaller in that band. The R — I 
colour provides an approximate way to distinguish between the two 
companion spectral types allowed by a given I band excess, while 
the J — H colour is not as discriminating (see dashed contours in 
Fig.H. 

In actuality we do not rely on the diagrams in Figs. [2] to de- 
termine the companions' spectral types and their limits. Rather we 
predict the primary's I (J) band absolute magnitude using the dis- 
tance (Table [TJ and, together with the observed one we derived the 



secondary's absolute I (J) magnitude, which we then convert into 
a spectral type using the table in Appendix [Cl In this way we do 
not rely on the T c g — My relation characteristic of a given cool- 
ing curve. However, the companion spectral types listed in Tables[8] 
and [9] are consistent with those that would be derived using Fig. [2] 
within one spectral subtype, except for the two objects with com- 
panions brighter than spectral type K0V, of which we speak more 
in §§[§3]and[12] 



6 RESULTS 

Those objects for which A(V — I) > (Ja{v-i) ( or » equivalently, 
A(V — J) > (Ja{v-j)) are considered as cases where a compan- 
ion is detected. We note that this difference is the same as what we 
call the / (J) band excess in Fig. [2] because the theoretical colours 
are normalised to the observed (dereddened) V magnitudes. Red- 
denings, intrinsic V — I (V — J) and R — I (J — H) colours, I (J) 
band excesses (A(V — I) or A(V — J)), companions' absolute I 
(J) band magnitudes and companions' spectral types are listed in 
Tables [8] l[9}. The / (J) band excesses are plotted as a function of 
stellar temperature in Fig.[3]l|4]l. 

The reddenings derived by comparing observed and (single 
star) predicted B — V colours are compared to those derived from 
nebular methods in Fig.|5] Most reddenings compare well. The five 
exceptions are A 57 and A 79, whose high derived stellar redden- 
ings are likely due to bright companions affecting the B and V 
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I Excess (Mag), R-l Color (Mag) 

M8M7M6 M5 M4 M3 M2 MO K5 KO G5G0 F5 



J Excess (Mag) J-H Color (Mag) 

M8M7M6 M5 M4 M5 M2 MO K5 KO G5 GO F5 





Companion Moss (M ) 



Companion Mass (M©) 



Figure 2: The / and J band excess (solid contours) and I — R and J — H colours (dotted contours) that are produced by main sequence companions (M9V 
to F5V) as a function of companion mass and primary's absolute V magnitude. These plots are obtained by assuming that the central stars are on the cooling- 
curve part of the HR diagram and cannot be used to determine the spectral type of the companion of objects which have not yet turned the "knee" on the 
post-AGB evolutionary track. The / and J band excess correspond to the A(V — /) and A(V — J) listed in Tables[8]and[9] 



Name 


E(B - V) 


(V - I)o 


{R-I)o 


A(V-I) 


M I2 




Comp. spec, type 


A7 


0.02 ± 0.02 


-0.35 ± 0.08 


-0.19 ±0.06 


-0.01 ± 0.08 


> 9.60 




Later than M4V 


A 16 


0.13 ±0.03 


-0.12 ±0.11 


0.02 ± 0.09 


0.22 ±0.11 


8.30 [9.24 -7.73 ] 


M3V[M3V-M1V] 


A 20 


0.07 ± 0.02 


-0.32 ± 0.08 


-0.16 ±0.06 


0.02 ± 0.08 


> 7.29 




Later than MOV 


A 28 


0.02 ± 0.07 


-0.35 ± 0.24 


-0.19 ±0.18 


-0.02 ± 0.24 


> 8.18 




Later than M2V 


A31 


0.00 ± 0.02 


-0.29 ± 0.08 


-0.14 ±0.06 


0.06 ± 0.08 


> 9.14 




Later than M3V 


A 57 


0.46 ± 0.04 


-0.03 ±0.15 


0.15 ±0.12 


0.29 ±0.15 


5.51 [6.56-4.79] 


K3V [ K7V - G8V ] 


A71 


0.38 ± 0.02 


-0.31 ±0.08 


-0.01 ± 0.07 


0.05 ± 0.08 


> 10.56 




Later than M5V 


A 72 


0.04 ± 0.04 


-0.35 ±0.14 


-0.15 ±0.14 


0.01 ± 0.14 


> 6.84 




Later than K8V 


A 79 


1.20 ±0.01 


-0.21 ±0.10 


0.42 ±0.10 


0.15 ±0.10 


3.12 [4.44-2.66] 


F6V [ G6V - F3V ] 


A 84 


0.11 ±0.03 


-0.22 ±0.10 


-0.08 ± 0.08 


0.12 ± 0.10 


10.00 [ 11.91 - 


9.28 ] 


M4V [ M6V - M3V ] 


DeHt5 


0.07 ± 0.03 


-0.22 ±0.11 


-0.08 ± 0.09 


0.10 ±0.11 


> 9.28 




Later than M3V 


EGB 1 


0.20 ± 0.01 


-0.29 ± 0.06 


-0.07 ± 0.05 


0.07 ± 0.06 


10.05 [ 11.96- 


8.91 ] 


M4V [ M6V - M3V ] 


EGB 6 


0.03 ± 0.01 


-0.34 ± 0.03 


-0.17 ±0.03 


0.01 ± 0.03 


> 9.79 




Later than M4V 


HaWe5 


0.14 ±0.03 


-0.25 ±0.10 


-0.08 ± 0.08 


0.06 ± 0.10 


> 11.21 




Later than M5V 


HDW3 


0.23 ± 0.01 


-0.32 ± 0.04 


-0.06 ± 0.05 


0.03 ± 0.04 


> 10.35 




Later than M4V 


HDW4 


0.04 ± 0.02 


-0.24 ± 0.07 


-0.11 ±0.06 


0.07 ± 0.07 


13.05 [ 16.87- 


12.30] 


M7V [ M8V - M6V ] 


IsWe 1 


0.19 ±0.03 


-0.35 ± 0.09 


-0.11 ±0.08 


0.01 ± 0.10 


> 9.77 




Later than M4V 


IsWe2 


0.32 ± 0.04 


-0.32 ±0.16 


-0.04 ±0.12 


0.04 ±0.16 


> 9.38 




Later than M4V 


JnEr 1 


0.00 ± 0.01 


-0.37 ± 0.06 


-0.21 ±0.05 


-0.01 ± 0.06 


> 10.86 




Later than M5V 


K 1-13 


0.00 ± 0.07 


-0.42 ± 0.24 


-0.25 ±0.19 


-0.09 ± 0.24 


> 8.85 




Later than M3V 


K2-2 


0.02 ± 0.02 


-0.31 ±0.09 


-0.17 ±0.07 


0.01 ± 0.09 


> 7.40 




Later than M1V 


NGC 3587 


0.00 ± 0.04 


-0.42 ±0.15 


-0.23 ±0.12 


-0.07 ±0.15 


> 8.96 




Later than M3V 


NGC 6720 


0.00 ± 0.03 


-0.29 ±0.11 


-0.16 ±0.08 


0.06 ±0.11 


> 7.99 




Later than M2V 


NGC 6853 


0.00 ± 0.03 


-0.32 ±0.11 


-0.16 ±0.08 


0.04 ±0.11 


> 8.42 




Later than M3V 


PuWe 1 


0.08 ± 0.02 


-0.34 ± 0.09 


-0.14 ±0.07 


0.01 ± 0.09 


> 10.46 




Later than M4V 


Sh 2-78 


0.31 ± 0.01 


-0.26 ± 0.06 


0.00 ± 0.06 


0.10 ±0.06 


10.22 [ 11.32- 


9.72] 


M4V [ M5V - M4V ] 


Sh 2-176 


0.26 ± 0.09 


-0.30 ± 0.32 


-0.03 ± 0.25 


0.05 ± 0.32 


> 8.91 




Later than M3 


Sh 2-176 1 


0.23 ± 0.02 


-0.25 ± 0.09 


-0.01 ± 0.07 


0.11 ±0.09 


10.68 [ 12.87 - 


9.95 ] 


M5V [ M6V - M4V ] 


Sh 2-188 


0.31 ± 0.01 


-0.30 ± 0.05 


-0.04 ± 0.04 


0.05 ± 0.05 


10.46 [ 14.33 - 


9.71 ] 


M4V [ M8V - M4V ] 


Ton 320 


0.00 ± 0.04 


-0.38 ±0.14 


-0.21 ±0.11 


-0.05 ±0.14 


> 9.24 




Later than M3V 


WeDe 1 


0.06 ± 0.01 


-0.34 ± 0.03 


-0.16 ±0.03 


0.01 ± 0.03 


> 11.36 




Later than M5V 



These measurements are obtained by excluding the observations taken in night 3 



Table 8. /-band excesses (A (V — I)), companion absolute /-band magnitudes (Af/a) and spectral types (or limits) of our targets 



bands, A 71, EGB 1 and Sh 2-78, for which we suspect the nebular 
reddenings have uncertain values. 

We have 9 detections or marginal detections in the / band: 
A 16 (at the 2-a level), A 57 (at the 1.9-ct level), A 79 (at the 
1.5-cr level), A 84 (at the 1. 2-a level), EGB 1 (at the 1.2-cr level), 
HDW 4 (at the l-a level, but this is a PN mimic), Sh 2-78 (at the 



1.7-cr level), Sh 2-176 (at the 1.2-cr level, only for the measure- 
ment excluding night 3) and Sh 2-188 (at t he 1-cr level). A 79 wa s 
already known to have a F0V companion dRodrfguez et al.ll200ll) . 
The spectral type determined with our method is cooler because the 
reddening was overestimated due to the strong contribution of the 
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Name 


E(B - V) 


(V 


- J)o 


(J-H) 


A(V - J) 


M.J2 




Comp. spec, type 


A7 


0.02 ± 0.02 


-0.66 


±0.11 


-0.07 ± 0.20 


0.15 ±0.11 


9.65 [ 11.07 -S 


.95 ] 


M5V [ M8V - M5V ] 


A31 


0.00 ± 0.02 


-0.41 


±0.07 


0.15 ±0.03 


0.40 ± 0.07 


8.26 [ 8.57-7 


98] 


M4V [ M4V - M4V ] 


A 72 


0.04 ± 0.04 


-0.68 


±0.17 




0.16 ±0.17 


> 6.57 




Later than Ml 


A 79 


1.20 ±0.01 


-0.74 


±0.06 


0.02 ± 0.07 


0.10 ±0.06 


4.05 [5.17-3.45] 


G7V[K5V-G1V] 


DeHt5 


0.07 ± 0.03 


-0.24 


±0.14 




0.54 ±0.14 


8.80 [9.30 -8 


37] 


M5V [ M5V - M4V ] 


EGB 1 


0.20 ± 0.01 


-0.65 


±0.11 




0.18 ±0.11 


9.43 [ 10.52 -i 


.80] 


M5V [ M6V - M5V ] 


EGB 6 


0.03 ± 0.01 


-0.54 


±0.09 


0.38 ±0.10 


0.28 ± 0.09 


9.15 [9.67-8 


76] 


M5V [ M5V - M5V ] 


K2-2 


0.02 ± 0.02 


-0.71 


±0.09 


-0.06 ± 0.08 


0.01 ± 0.09 


> 7.84 




Later than M4V 


NGC 3587 


0.00 ± 0.04 


-0.93 


±0.18 




-0.12 ± 0.18 


> 10.12 




Later than M6V 


NGC 6720 


0.00 ± 0.03 


-0.63 


±0.22 




0.19 ±0.22 


> 7.97 




Later than M4V 


NGC 6853 


0.00 ± 0.03 


-0.66 


±0.11 


0.05 ± 0.06 


0.17 ±0.11 


8.81 [ 10.05-8 


.12] 


M5V [ M6V - M4V ] 


Ton 320 


0.00 ± 0.04 


-0.87 


±0.20 




-0.08 ± 0.20 


>9.92 




Later than M6V 



Table 9. J-band excesses (A (V — J)), companion absolute J-band magnitudes {Mj%} and companion spectral types (or limits) of our targets. All detections 
and limits are consistent with the results of the /-band excess (Table[8j 




-0.5 t i i , i I i I I i I I i i i i I , i i i □ 

50 100 150 200 

T eff (kK) 

Figure 3: The observed (dereddened) V — I colours of our targets (symbols; see Table[8j compared with the predicted V — I colours of single stars as a 
function of effective temperature (for a log g = 7.0 - solid line) 



companion in the V band (see §[5). We have therefore 8 detections 
out of 27 bona fide objects. 

In the J band we detected A 7 (at the 1.4-cr level), A 31 (at the 
4.7-0- level), A 79 (at the 1.7-ct level), DeHt 5 (at the 3.8-cr level, 
but this is a PN mimic), EGB 1 (at the 1.6-ct level), EGB 6 (at the 
3.1-0- level) and NGC 6853 (at the 1.5-er level). We have therefore 
6 detections out of 1 1 bona fide central stars with data. 

The spectral types of the companions implied by the detected 
excesses as well as limits of the non-detections are also listed in 
Tables [8] and [9] When a companion was not detected, we summed 
the I (J) band excess with the upper error bars to create an upper 
limit to the flux excess and determined an upper limit for the com- 
panion mass/spectral type in this way. We note that all limits and 
detections are fully consistent across the two detection techniques. 

The / band binary fraction is determined by the ratio of 8 ob- 
jects over 27 bona fide central stars (see §|2]for an explanation of 



what we mean by bona fide), or 30%. Only companions brighter 
than the spectral type M3-4V can be detected by our survey, where 
this limit was estimated by taking the median of the limits in Ta- 
bled] Within the J band group, the detection rate is 6 out of 11, or 
54%, in line with the study of iFrew & Parked d2007l . see also Frew 
2008). A look at Table [9] shows that most detections in the J band 
were fainter, in line with the expectation that the J band is a more 
sensitive method. Based on the detection and limits of Table [9] we 
guesstimate that this limit is M5-6V. The binary fraction for the en- 
tire sample, obtained from either the I or J band excess methods is 
44% (12 out of 27 bona fide PN), to a limit intermediate between 
the two methods. 

Before comparing our results with predictions, we would like 
to remark on the asymmetry of the distribution of V — I colours of 
our sample about the single star prediction line (Fig.|3j- Even elim- 
inating from that plot the 9 detections, which are all above the pre- 
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Figure 4: The observed (dereddened) V — J colours of our targets (symbols; see Table[9} compared with the predicted V — J colours of stars as a function 
of effective temperature (for a logg = 7.0) 




0.2 0.5 0.8 1.1 
nebular— derived E(B— V) 



Figure 5: A comparison of the reddening values obtained via a nebular 
method (Table[T}, and those obtained via our stellar method (Table[8j 



diction, the non-detections (21 objects) still preferentially lie above 
the single star prediction line (15 vs. 6 objects). We speculate that 
this points to the detection of a binary signal which is higher than 
the fraction we determined from individual objects above. To de- 
termine whether this is significant we would need a Monte Carlo 
simulation. One could counter argue that this asymmetry is caused 
by an incorrect placement of the prediction line, or, in other words, 
that the predicted colours are too blue. To shift the prediction line 
so that the non detected data points scatter symmetrically about it, 
would require a systematic shift of the synthetic V — I colours 
of about ~5% which is excluded by the comparison of synthetic 



spectra wit h data carried out o ver the years (Rauch, private com- 
munication; [Rauch!etaj]|2007l). 

Finally, we have demonstrated that the / band excess method 
necessitates extraordinary accuracy, such as that achieved for the 
current dataset, but is not usually encountered in unvetted data 
from the literature. For example, in Fig. [6] we compare the J-band 
excesses obtained by using the B, V and 7 data compilation of 
Bilfkova et with our dataset. Using all of their data and 

determining the stellar temperatures in the same way we have done 
for our sample, we see how the scatter around the single star pre- 
diction is dramatic. While their compilation was not aimed at de- 
tecting binarity, it can easily be seen that data from the literature is 
not generally suitable for this type of work. 



6.1 Photometric Variability 

Periodic variabilit y in central stars of P N usually denotes a short 
period binary (e.g jMiszalski et alj2009h . Variability should not in- 
terfere with the detection of an IR excess because the exposures in 
each filter are taken in a short seq uence and the periods o f these ob- 
jects are of the order of 12 hours jPe Marco et alj2009l) . Although 
it is expected that the colour of the short period bi nary should vary 
durin g the orbit, this is only marginally so (e.g., |Pe Marco et al.1 
2008). Variability can however increase the error on the absolute 
photometry because we average observations taken at more than 
one epoch. We have therefore monitored all targets for variability 
during photometric as well as non-photometric nights. The differ- 
ential magnitudes to 5 field stars are then plotted as a function of 
time. The central stars for which the scatter of the average differen- 
tial magnitudes in each filter is larger than the mean errors on these 
magnitudes are labelled as possible variables. Further, we retain as 
likely variables those of these stars for which the variability trends 
are the same in all filters. We will continue monitoring these stars 
to determine if a period is present. The variable stars are: A 57, 
A 72, A 84, IsWe 2, Sh 2-78 and Sh 2-176. Four of the 6 variables 
have an / band excess, as one might expect if the variability de- 
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Fi gure 6: Left panel: the same figure as Fig. [3] except here we compare our results (with no error bars, for clarity) with results obtained using the data compiled 
bv lBilfkova et al.l (2012, right panel), showing how unvetted data has much larger random errors 



notes binarity. Further comments on the variables are reported in 



7 COMPARISON OF THE OVERALL PN BINARY 

FRACTION WITH THE OVERALL MAIN SEQUENCE 
BINARY FRACTION 

We recognise that the binary fraction determined here is prelim- 
inary, because of the small sample size. However, here we com- 
pare it with the prediction from the current PN evolution scenario, 
whereby PN derive from all ~1-8-Mq stars, whether they are sin- 
gle stars or are in binaries. As already explained, it is this prediction 
that we are testing. We expect that if PN derive from a binary in- 
teraction more often than for the current evolutionary scenario, we 
would find a larger fraction of short- and intermediate-period bina- 
ries in the central star population. This leads to the prediction that 
the overall binary fraction, i.e., where we consider binaries at any 
period, should also be larger in the binary hypothesis than in the 
standard scenario. 

The progenitor population of the PN from which we have 
drawn our targets is the main sequence stars with masses be- 
tween ~0.9 and 8 Mq. We can then use t he ma in sequence pop- 
ulation binary fraction of Raghavan et al. (2010) to represent the 
binary fraction of the progenitor population of our PN. We use 
(50±4)%, which is the percentage of double and multiple star sys- 
tems with primaries in the spectral type range F6V-G2V (masses 
in the range 1.27 - 1.03 Mq; see Table ICTV This is reasonable 
on the g rounds that the median p rogenitor mass of today's PN is 
1.2 M dMoe & De Marcoll2006t). This frac tion is lower than that 
determined by Duqu ennov & Mayor] dl99 ll 57%) who accounted 
for a larger incompletion bias, but is likely to be more accurate. 
This fraction includes any main sequence companion down to the 
planetary regime at any separation. 

7.1 Accounting for completion effects 

Before we can carry out a comparison we need to increase the PN 
binary fraction to (i) include the companions that have not been 
detected because they are too faint; (ii) include wide companions, 
because our observed sample does not contain resolved binaries by 
design. Finally, we need to (iii) account for the fact that some main 



sequence close binaries go through a common envelope interaction 
on the RGB. Those systems will become short period binaries with 
very lo w mass envelopes (the prim ary will become a subdwarf O or 
B star; iMorales-Rueda et "alll2003h and are unlikely t o ever ascend 
the A GB on the grounds of low envelope masses dDorman et al.l 
1993); if they do, they will suffer a second common envelope al- 
most immediately which would result in a sub-luminous object, 
highly unlikely to ever make a visible PN. 

(i) To account for companions with an M3-4V spectral type or 
fainter (mass ^0.33-0.24 Mq), we use Table 18 of lRaghavan et al.l 
d2010h that lists all the companion spectral types of the detected 
binaries. We determine that the fraction of main sequence binaries 
with a companion spectral type of M3-4V or later is 41-31% (73- 
56 of 179 main sequence binaries for which the companion spec- 
tral type is known, where we count as undetectable also 4 hot and 
evolved companions). We therefore increase our /-band derived PN 
binary fraction (30%) by a factor of 1 .69- 1 .45 to account for the un- 
detected faint companions. 

To account instead for companions with brightness equal to, or 
fainter than M5-6V (the limit we estimated for the J band dataset) 
we would have to multiply the J-band-derived PN binary fraction 
(54%) by a factor of 1.28-1.19. 

(ii) The ground-based spatial resolution of our observations is 
probably only slightly smaller than the median seeing of 1.3 arc- 
sec corresponding to a projected separation of ~1040 AU for the 
median distance of our sample of 0.80 pc. To de-project this sep- 
aration we would have to divide the projected separation by a fac- 
tor of 0.8, which accounts for the systems' random phase and ran- 
dom orientation. However, if we also account for an eccentricity 
distribution similar to the main sequence's, this factor approaches 
unity, because systems spend more time near apastron in eccentric 
systems. Not knowing the eccentricity distribution of central star 
binaries but presuming a degree of circularisation to have taken 
place, we de-project the projected separation using a range of fac- 
tors (0.8-1.0). The median separation is then converted to a period 
in the range ~13000-18 000 years (log(P/days)=7. 11-7.25), using 
Kepler's third law and a total system mass of 0.9 Mq. This pe- 
riod limit is then decreased to log(P/days) = 6.71 — 6.85, to 
account for the fact that, due to mass-loss and angular momentum 
conservation, the orbit of a typical binary will widen by a factor 
of 2-3 (we used 2.5) and so some of the main sequence binaries of 
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Ragh avan et all J2010h will become resolved central star binaries in 
our sample because of orbit widening. 

A similar argument can be applied to the J-band sample of 
1 1 bona fide PN, which have a median distance of 0.65 kpc, with 
a resulting median projected separation of 845 AU, resulting in a 
period corrected for orbit widening of log(P/days) = 6.57 — 6.71. 

(iii) Finally, the main sequence binary population with pe- 
riod shorter than log(P/days)=2.43 never ascend the AGB. This 
value was calculat ed using a radius on th e RGB of 100 R0(see, 
e.g., see figure 5 of |Pe Marco et all l l201ll) ). and a maximum tidal 
capture radius of 2 stellar radii (using a total system mass of 
1.5 Mq). T he maximum tidal captur e rad ius for the RGB was ob - 
tained from I Villaver & Livid J2009h and iMadappatt et all j2oTTb . 
At this point we integrate under the normalised m ain sequence bi- 
nary period distribution of Rag havan et afl feOld their Fig. 13) us- 
ing log(P/days)=2.43 - 6.78 as limits (or 2.43 - 6.65 for the J-band 
sample; where we have taken the average of the period separation 
ranges to determine the upper limits), obtaining a fraction of 0.65 
of the total (or 0.64, for the J band sample): these are the systems 
we detect. This translates into a second factor of 1.54 (or 1.56), 
which we multiply by the central star binary fraction so as to in- 
clude binaries at all separations and thus make it comparable to the 
main sequence one. 

(iv) The only bias that is impossible to account for, is that due 
to central stars whose companions are bright enough to contribute 
to the V and even B bands. As we have explained in §|5] these will 
result in reddenings that are artificially large, because the bright 
companions tend to make their system's colours redder. Once this 
over-large reddening is applied to all the bands, the binary SED 
is rendered bluer than it should be and the I or J band excess is 
necessarily reduced. The mismatch between reddening curve and 
SED of the companion can decrease the red/IR flux excess below 
detectability or reduce it such that it is impossible to match it to a 
companion spectral type. Interestingly, only 6 objects have stellar- 
derived reddening values that are larger than the nebular-derived 
values, the most noticeable being A 79, a known binary and A 57 
which has a large I band excess. For the other 4 objects (A 7, A 16, 
A 28 and Sh 2-188), these discrepancies are however always within 
the uncertainties. We therefore conclude that either the reddening 
comparison is not reliable due to high errors in nebular-derived red- 
denings, or that our sample tends not to include companions much 
hotter than spectral type K0V. This makes sense in view of the fact 
that only a minority of companions are expected to be that hot). 
We leave this bias unquantified in the knowledge that it would be 
at most a few percent. A spectroscopic follow-up will resolve this 
issue. 

7.2 The debiased PN binary fraction and its uncertainties 

The de-biased central star binary fraction obtained through the / 
band photometry of our sample is 67-78%. Using the J band data 
alone we obtain a fraction of 100-107%. The ranges account for 
the uncertainty on the factor to account for the unobservable faint 
companions. The error bars on these estimates are difficult to esti- 
mate for the moment and we defer this exercise to the next paper 
in this series, which will increase the sample. However, if all limits 
were accurately accounted for then the I- band and J-band frac- 
tions should be the same. We therefore tentatively estimate a 10- 
30% error bar. We preliminarily conclude that the PN binary frac- 
tion is higher than the main sequence binary fraction of (50±4)%. 

We can also estimate the PN binary fraction after eliminat- 
ing from each of the I and J samples the one detection with the 



smallest statistical relevance. By doing so, we would obtain / and 
J binary fractions of 58-68% and 84-90% respectively. 



7.3 Comparison of the short-period PN binary fraction with 
the main sequence binary fraction 

The fraction of PN that surround short-period, post-comm on enve- 
lope binaries is 15-20% dBondl200d : lMiszalski et alj2009h . This is 
a very large fraction when we consider that only a very small mi- 
nority of main sequence binaries would suffer a common envelope 
on the AGB. 

The fraction of main sequence binaries that go through a com- 
mon envelope on the AGB, resulting in post-common envelope cen- 
tral stars of PN should be quite small: only those companions that 
escape capture on the RGB, but then are successfully captured on 
the AGB will become post-common envelope central stars. These 
are companions that, while on the main sequence, have periods in 
the range log(P/days) > 2.7 — 2.8 and therefore account for only 
~ 1 % of all main sequence binaries. The lower period limit was cal- 
culated in § [7j while the higher period limit was obtained by con- 
sidering that for a successful AGB capture the orbital separation 
has t o be smaller than approximately twice the maximum AGB ra- 
dius JVillaver & Liviol 2009; Nor dhaus et alj20ld : lMadappatt et all 
l201lh . The maximum AGB stellar radius for a 1 .2-Mq star was 
calculated to be ~300 R© foe Marco etai]|20 111 with the usual 
adjustment for the orbital widening due to mass loss, see §|7). Even 
exaggerating the maximum AGB radius (600 Rq) and the maxi- 
mum capture radius to 5 stellar radii, we get log(P/day) = 3.9 and 
a fraction of 5%. We would therefore predict a similar fraction of 
post-common envelope central stars, contrary to the observations. 



7.4 Comparison of the PN binary fraction with the white 
dwarf binary fraction 

We note, finally, that the PN binary fraction, even before detec- 
tion biases are accounted for, appears to be higher than the white 
dwarf binary fraction of 25% dHolberg 2009, the WD binary frac- 
tion should be a few points higher, if we could readily detect white 
dwarfs around bright main sequence stars - Sirius-like systems). 
One reason could be that the white dwarfs derive from a population 
that includes lower mass stars, those that do not develop a PN due 
to long transition times. Such low er mass population would natu- 
rally have a lower binary fraction ( Rag havan et al] |201Ch . Another 
reason could be that PN form preferentially around binaries, the 
hypothesis we are trying to test. To discern between the two expla- 
nations of the discrepancy between the white dwarf and PN binary 
fractions we need better constraints on white dwarfs and central 
stars masses dLiu et alJl995l : lLiebert et alj2005klGesicki & Ziilstral 
120071) . 



8 NOTES ON REMARKABLE INDIVIDUAL OBJECTS 
8.1 Abell 7 

This star was detected in the J band (at the 1.4-cr level) and the 
implied companion would be an M5V star. Companions dimmer 
than M4V are excluded by the I band photometry. The central star 
has a red-dwarf companion 0.91 arcseconds away from the central 
star (Ciardullo et al. 1999) resolved by the Hubble Space Telescope 
(HST), which warrants further investigation in order to determine a 
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photometric parallax; only an upper limit of K2V could be deter- 
mined from the HST data of Ciardullo et al. (1999), while Frew 
(2008) determined a spectral type of M4V, similar to our own. 

8.2 Abell 16 

We detected an M3V companion in the I band at the 2-a level. 

8.3 Abell 31 

Ciardullo et al. (1999) give an upper limit to the distance of 440 pc 
(based on photometry of the resolved companion), mildly inconsis- 
tent within the uncertaint ies with a direct tr igonometric determina- 
tion of D = 56819Q 1 pc dHarris et al1l2007h . Our I and J photom- 
etry su£gj«tjin_M4y_£ompanion, in agreement with Frew (2008) 
and lCiardullo et alj {l999), who used the HST to determine a sep- 
aration of 0.26 arcsec for a companion with spectral type later than 
M4V. This companion may never have interacted with this PN. The 
morphology of this PN is likely round, in agreement with the lack 
of interaction. 



8.4 Abell 57 

The /-band excess detected for this object exceeds the uncertainty 
by a factor of two and the spectral type of the companion is K3V. 
Such a bright companion may have contributed in the V band. In 
fact, the reddening derived from our method is higher than the red- 
dening derived from nebular observations, (cf. Tables[T]and[6}. We 
also note that we used the limiting temperature of 60kK for this star. 
A higher temperature would lower the / band excess slight ly, with 
a temp erature of 150kK reducing it to 0.25 mag. iMiszalski et all 
d20 1 lh discovered this t o be another example of a EGB 6-like cen- 
tral star lFrew & Parker! d2010l) . implying the presence of a compact 
disk either around the central star or its companion, lending strong 
support to this star being a physical binary. This central star appears 
to be variable at the 0. 1-mag level from our differential photometry, 
while our absolute photometry reveals at most a variability of a 2% 
based on two epochs of data. 

8.5 Abell 72 

As was the case for A 57, the effective temperature of the central 
star of A 72 can be constrained to be larger than 100 kK by the 
Zanstra method. We therefore adopt a lower limit of 100 kK due 
to the fact that our own UVES spectrum suggests a temperature in 
excess 120 kK. There is effectively no difference in the calculated 
excess for temperatures higher than 100 kK. From I and J band 
photometry we can state that the companion, if present, would be 
fainter than the Ml V spectral type. We suspect this star to be vari- 
able at the 0.05-0. 1-mag level. 

8.6 Abell 79 

This star is known to have a cool spectral type dRodrfguez et al.1 
l200lh possibly denoting that it is a binary since the ionising source 
of the PN would otherwise be missing. When the cool companion 
dominates the light of the system, we can assume that it contam- 
inates the V or even the B bands. As a result, the reddening will 
be too large. We determined E(B — V) = 1.2 while from near-IR 
data a value of 0.62 is derived (Frew 2008). Using our reddening, 



all the photometry is de-reddened excessively, with resulting de- 
rived V — I and V — J colours that are too blue and a reduced / 
and J band excess. That is why the determined companion spec- 
tral type (F6V from / band photometry and G7V from J band 
measurements) are redder/fa inter than the one that w as determined 
from spectroscopy (GOV) by Rodrigue z et alj J200lh . Note how the 
result from J-band photometry is closer to the spectroscopically- 
derived spectral type as predicted in §|5] In addition, we cannot use 
the contour plots in Figs.[3]and|4]because we do not have the central 
star's My value. 

The / band excess for this star is large (0.15 mag) but it is 
affected by a considerable uncertainty. The reason is that the / band 
magnitude of this star is 0.2 mag brighter in night 4 and in night 3. 
While we could exclude the N3 measurement on the grounds that 
the night was only partly photometric, our logs show that that part 
of the night was completely clear. In addition, the magnitudes in 
all other bands are very similar across the two nights. We therefore 
keep both measurements and accept the large uncertainty. 



8.7 Abell 84 

The central star of A 84 has a companion detected in the J band 
at the 1.2-ct level, making it a marginal detection. The differential 
photometry shows a clear dimming trend over the 8 nights of data 
at the 0. 1 mag level. 



8.8 DeHt 5 

This larg e emission ne b ula (a lso known as DHW5) was first iden- 
tified by iDengel et al.1 dl980h as a likely PN. It has an unusual 
irregular morphology, and shows quite a different appearance in 
[O III] li ght, compared to its morphology in [N II] and [S II] 
emission (Rosado & Morendfl99ll : lTweedv &~Kw itter 1996]). Pre- 
vious researchers hav e assumed that DeHt 5 is a bona fide PN. 
iBannister et alj d2003h had previously noted the proximity of its 
central star, WD 2218- 706, to the giant molecular cloud complex 
described by Kun Jl998l) . and stated that the star may lie in an area 
where the ISM is dense. However, these authors did not consider 
the possibility that the optical nebula might be in fact ionised ISM. 
Tweedy & Kwitter (1996) also noted the extensive diffuse emission 
around the bright nebula, suggesting this might be ionised interstel- 
lar material, but assumed that the core was a true PN. Frew (2008) 
found that the nebular and stellar velocities are different, suggesting 
that the n ebula is unrelated to its ionising star, a situation similar to 
PHL 932 dFrew et al1l2010l) . The overall body of evidence is again 
in favour of the ionised-ISM interpretation, rather than a PN. The 
morphology is not typical of an evolved PN, and the gas is consis- 
tent with being ionised ambient material, as the systemic velocity 
agrees with the CO velocities of widespread neutral gas at ~360 
pc, the distance of the star as determined by Benedict et al. (2009). 
The narrow line width of the nebular gas is also consistent with the 
ionised ISM interpretation. Furthermore, the evolutionary position 
of the star in the HR diagram is not in agreement with a post-AGB 
track that is consistent with the timescale of PN evolution. All of 
the evidence is consistent with DeHt 5 being a HII region ionised 
by WD 2218+706. The central star has a detected J-band excess, 
consistent with an M3V companion, but we flag that this is one of 
two objects with one epoch of measurements only. 
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8.9 EGB 1 

This somewhat irregular, one-sided nebulosity is traditionally clas- 
sified as a PN. It is poorly studied, but spectra show moderately 
bright [O III] emission. However, the lack of a limb-brightened 
bow shock (expected for such a on e-sided nebula wit h an off-centre 
central star), the small line-width dLopez et alj|2012f) . the systemic 
nebular velocity close to zero in the local standard of rest (LSR), 
and the moderate excitation, lead to questions regarding its status. 
It may be another case of ionised ISM, but more observations need 
to be obtained before a definitive answer can be given and there- 
fore we have kept it as a PN for the time being. Additionally, the 
ionising star is a periodic photometric variable (T. Hillwig, private 
communication), consistent with binarity. We detected an / and J 
band excess at the 1.2 and 1.6-<r levels, respectively, consistent with 
an M4-5V companion. The hypothesis that the companion to this 
central star could be a white dwarf is inconsistent with the R — I 
colour. 



8.10 EGB 6 



iLiebert et alj dl989h found an unresolved, dense emission neb- 
ula surroun ding the central star, w ith a tiny ionised mass of 
~10 -9 M m . IZuckerman et alj < fl99lt) found near-IR evidence for a 
probab le red dwarf companion, confirmed by iFulbright & Lieber3 
dl993l) via JHK photometry o f the central star. Th e companion 
has been imaged with HST (see lBond et all 19931 and lBondll 19941) 
at a projected separation of 0.18". Furthermore, the unresolved 
emission nebula correspond s in position to the red dwarf and not 
the central star jBondll 19941) . suggesting there may be a small disk 
of material around the cool companion (possibly accreted from the 
PN), and photoionised by the central star. This is the archetype of a 
class of PN centr al stars with compact, unresolved ionised nebulae 
(or disks), which Frew & Parked (2010) refer to as EGB 6-like cen- 
tral stars. We detected a 3-cr J band exc ess consiste nt with an M5V 
companion, in line with the results of iBond et alj dl993l) . but we 
note that the J — H colour is mildly inconsistent with this spectral 
type, as it appears too red. 



8.11 HDW 4 



Also d e signated HaWe 6, t his faint nebula was noted bv lHartl et al.1 
dl983l) . lHarriset al. (2007) determined a trigonometric distance to 
the ionising star of only 2091 is pc. Napiwotzki (1999) performed 
a NLTE model atmosphere analysis, determining T e g = 47,300 
± 1700 K and logg = 7.93 ± 0.16, and estimating a distance 
of 250 pc. The relatively low temperature and high gravity of the 
white dwarf suggests a long cooling age far greater than any f ea- 
sible PN lifetime. Referring to Table 1 of Berge ron et al.1 dl995h . a 
cooling age of 3-4 xlO 6 years is suggested. The absolute magni- 
tude is My = 9.43lo^g which is considerably fainter than any PN 
nucleus. The very low ionised mass of the nebula (5 x 10~ 3 M B ; 
Napiw otzki 1999) also rules out a PN interpretation. iNapiwotzlg 
( 1999) speculates that the nebula might be a shell produced by an 
ancient nova outburst (the mass is approximately right), but the star 
shows no sign of any cataclysmic variable features in its spectrum 
dNapiwotzki & Schon berneg|l995l) . Napiwotzki (1999) obtained a 
spectrum of the HII region around this star and deduced an upper 
limit to the expansion of 2u cxp < 47 km s - . Hence an expanding 
nova shell hypothesis is ruled out. The expansion limit is consis- 
tent with ambient ISM, so the nebula may be another case of ISM 
ionisation by the unrelated hot white dwarf WD 0533+555 (Frew 



2008). We report a very marginal l-cr detection in the I band, con- 
sistent with an M7V companion. However, this object was observed 
only once and that observation may have been affected by thin cir- 
rus (TableO. 



8.12 Ha We 5 

For the 'central' star, Napiwotzki (1999) determined T cB = 38,100 
± 1500 K and log g = 7.58 ± 0.20 from an NLTE model atmo- 
sphere analysis. The estimated distance is ~420 pc. This small faint 
nebula shares characteristics in common with HDW 4, where the 
large evolutionary age of the star is greatly at odds with the exis- 
tence of a remnant PN. iNapiwotzkH d 19991) also suggests that this 
may be an old nova shell and estimated an ionised mass of only 
2x 10 -4 Mq. A red Digitized Sky Survey image shows a vaguely 
PN-like elliptical nebula surrounding the ionising star. The small 
size and faintness of the object would be remarkable if it was a real 
PN at the nominal distance. Instead, it is more likely to be be a 
wisp of ionised ambient material. We note that the second epoch of 
this object was suspected non-photometric. The fluxes in the four 
bands for that exposure were indeed fainter than the other two by 
0.03, 0.23 0. 15 and 0.01 mag. While clouds are grey and one would 
expect a similar magnitude drop in all bands, it is possible that the 
high cirrus cloud may have changed over the stretch of time over 
which we took the 4 exposures. We therefore did not consider the 
exposures in this epoch for absolute photometry. On the other hand, 
the other two epochs differed by as much as 0.02 and this may be 
an indication of intrinsic variability. Our lack of I band detection 
imposes a limiting companion spectral type of M5 V. 



8.13 IsWe 2 

We are reasonably confident that this star is variable at the 3-4 per- 
cent level, although the clear trend in the B, V and R bands is not 
clearly evident in the I band (Table IaTT i. However, we noted that 
the I band photometry taken in the second epoch was 30 minutes 
later than the other filters, which may have contributed to the lack 
of a similar trend. From / band photometry we impose a limit of 
M4V on a companion's spectral type. The error on the average of 5 
epochs is possibly larger than it should be because of the variability. 



8.14 K 1-13 

This is a faint elliptical / bipo lar nebula also known as Abell 25. The 
ISchlaflv & Finkbeineij d201 lh asymptotic reddening in this direc- 
tion is E(B — V) = 0.03, so the central star is essentially unred- 
dened. We adopt an int egrated Ha flux o f logF = —11.95 from 
iFrew et alj d2012t) . and iKaler et ail (l990) gave an upper limit for 
the He II F(A4686) flux of 0.5 x FQIP). Using V = 18.42 for the 
central star, we obtain a hydrogen Zanstra temperature of 49 kK, 
which is a lower limit. The He II flux limit can be used to de- 
rive an upper limit on the temperature of 92 kK for this optically 
thick PN. We have t herefo re adopted a temperature of 80±30 kK. 
We note that I Abelll Jl966l) gave V = 18.94, B — V = -0.14, 
U — B = —1.15, so his l/-band magnitude is quite different to 
ours. We impose a limiting spectral type of M3 V from I band pho- 
tometry. 
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8.15 K 2-2 



The ionising star has a rather low mass dNapiwotzkill 19991) and the 
inferred post-AGB lifetime is not consistent with the kinematic age 
of the nebula. This object may turn out to be another case of ionised 
ISM. Like EGB 1, the nebula has a very low line width an d a sys- 
temic nebular velocity close to zero LSR dLopez et al.l2012l) . Alter- 
natively, this P N may be a senile ex ample of a common-envelope 
ejection event. [Afsar & BoncJ d2005l) have suggested that the cen- 
tral star is a likely close binary, based on radial velocity variability. 
Our J band photometry imposes a limiting spectral type for a com- 
panion of M4V. 



8.16 NGC3587 

This nebula is also k nown as M 97. The /-band magnitude of 
ICiardullo~et al. ( 1999) is fainter than ours by about one magnitude. 
We suspect it to be erroneous based on the fact that all available 
photometry (Sloan Digital Sky Survey, 2MASS and Spitzer/IRAC) 
is in line with our B, V and I values. We derive a companion spec- 
tral type limit of M6V from J band photometry. 



8.17 NGC6853 

This beautif ul, bright PN, also kn o wn as M 27, is at a di stance of 
405 pc l lBenedict et alj|2009l) . Iciardullo et all dl999l) discuss 
the companion to the central star discovered by |Cudwortrjdl973l) . 
Assuming a physical association, a tentative spectroscopic parallax 
of 430 pc is derived, in agreement with the trigonometric distance, 
but a radial velocity measurement for the companion is needed. We 
detect an M5V companion with J band photometry at the 1.5-er 
level. 



8.18 PuWe 1 

The effective temperature of the central star of PuWe 1 was a 
weighted average of the determination by Napiwotzki (1999), Good 
et al. (2004) and Gianninas et al. (201 0). Harris et al. ( 2007) give 
a trigonometric distance of 365^^ pc. Iciardullo et al] dl999h dis- 
cuss the nature of a pair of stars (likely a binary) that are located 
at 5.2 arcsec from PuWe l's central star. This pair may be a wide 
companion to the central star, in which case this system would be 
a triple. However, the proper motion of the pair appears different 
to that of the central star, which would make the central star not 
associated with it. Our photometry excludes the optical pair. Our 
/ band imposes a limit for a companion spectral type of later than 
M4V. 



8.19 Sh 2-78 

This star exhibits variability at the 0.3-mag level over the course of 
the 8 nights. The photometric magnitudes reveal only a ~0.05-mag 
brightening between nights 4 and 5. We detected an M4V compan- 
ion in the / band at the 1.7-cr level. 



8.20 Sh 2-176 

This object was detected to be variable at the 5% level from 
our monitoring dataset. Our three epochs of photometric observa- 
tions reveal that indeed the object brightened between the first two 
epochs and dimmed thereafter. Although the data in all four filters 



tell a consistent story, the third epoch data are affected by a larger 
uncertainty and the B band magnitude is quite a bit dimmer than 
the other filters. On the grounds of variability and larger errors we 
excluded the third epoch from the calculation of the average mag- 
nitudes (Table [6). By doing so the / band excess is increased and 
the error is greatly decreased, resulting in a 1.2-ct detection of an 
M5V companion (see the two entries in Table[8]l. 

8.21 Ton 320 

The effective temperature of this central star of Ton 320 was deter- 
mined as the we ighted mean of the estimates of Good et al. (2004), 
iTremblav et alj J201 1 ) and Gianninas et al. (2010), where the grav- 
ity was obtained from Gianninas et al. (2010). The ionising star 

of this very faint PN (also known as TK 1) has a 2MASS K-band 

1 — 11 -1 

excess (Holberg & Magargal 2005) but this may be spurious since 

the error bars on that value are large. The most stringent limit on 

the companion spectral type is imposed by our J band observations 

(M6V). 



8.22 WeDe 1 

There is a faint wide companio n to the central star of the PN 
WeDe 1 (also known as WDHS h lWeinberger et aljl983h . though 
it is likely to be an optical companion. Our / band photometry can- 
not detect a companion fainter than spectral type M5V. 



9 CONCLUSIONS AND DISCUSSION 

In this work we have started a survey of the 2-kpc, volume lim- 
ited sample to determine the binary fraction via a technique able 
to detect / and J band flux excess due to the presence of a cool 
companion. For this first survey we have selected 30 central stars 
of PN, of which we later determined three to be PN mimics. Of the 
remaining 27 central stars we have detected 8 with a possible I- 
band excesses (of which one was a previously known binary). For 
12 of these objects we collected the best J band photometry from 
the literature and determined that 6 out of 1 1 bona fide central stars 
have a flux excess. In total we have detected an excess in 1 1 of the 
27 central stars in either of the two bands. 

For the I band survey we calculate a detected fraction of 30%, 
which when de-biased to account for undetected objects, to result 
in an unbiased binary fraction of 67-78%. There are three principal 
de-biasing factors. The first to account for the fact that we do not 
detect companions fainter than M3-4V. The second, to account for 
the fact that we do not detect wide binaries by design. The third, to 
account for the fact that main sequence binaries that go through a 
common envelope on the RGB do not ascend the AGB and do not 
become PN. These factors have to be included in order to compare 
the PN binary fraction with the main sequence one. 

From the J band survey we calculate a binary fraction of 54% 
which, when debiased to account for binaries at any period and all 
companions, becomes 100-107%, clearly showing that the error bar 
due to the low number statistics is high. The discrepancy between 
the I and J band fractions is likely an effect of the low number 
statistics. Since the binary fraction determined using / band pho- 
tometry relies on more data than using the J band photometry, it 
is likely that the I band fraction is more accurate (even if the / 
band is less sensitive ). However the J band fraction agrees with the 
preliminary work of iFrew & Parked d2007h . who used 34 objects 
(although their detection limits were poorly quantified). 



The Binary Fraction of Planetary Nebula Central Stars I. 



17 



Thus debiased we can compare these fractions with 
that of F6V-G2V m ain sequence star binaries of (50±4)% 
( Ragh avanetaT1l20iCh . We preliminarily conclude that there may 
be an overabundance of central star binaries, compared to the puta- 
tive parent population. 

We have also noticed how, starting from the main sequence 
binary fraction and period distribution we would expect only a few 
percent of central stars of PN to be post-common envelope binaries, 
whereas previous surveys detected 15-20%. This discrepancy can 
only be reconciled within the binary hypothesis. 

Here we have demonstrated that the binary fraction is an elu- 
sive quantity. Even the best, most accurate observations cannot eas- 
ily detect faint companions and the only way to reduce the error is 
by analysing a much larger sample with similarly accurate photom- 
etry. Such high accuracy observations are extremely hard to obtain 
because of the need of photometric weather and the importance to 
obtain several observations on different nights. In the second in- 
stalment of this paper we will use the same technique with an addi- 
tional dataset and combine the results. Our goal is to analyse with 
similar accuracy the entire 2.5-kpc volume-limited sample of Frew 
(2008), comprising approximately 250 central stars. 

Finally, we compare the binary fraction determined here with 
two alternative prediction s. The first one wa s done on the basis 
of PN mor phology alon e. |Pe Marco & Sokerl d201lf) . revising the 
scheme of Soke predicted that about 60% of all central 
stars have interacted with a stellar companion. Here we measured 
that 48% of central stars have a companion closer than 1300 AU 
(30% x 1.45-1.69, where we took the middle of the range; § 17.11 1. 
These two numbers are easily the same within the uncertainty, but 
only if all of the companions we detected (and those we did not 
detect but accounted for) are closer than ~100 AU, i.e., have inter- 
acted with the central star. The comparison is possibly closer using 
the J band sample, whereby 67% of all central stars are binaries 
with any separation to a distance of ^1000 AU. 

Another prediction was that of iMoe & De Marcol ( feOHh . who 
used a population synthesis analysis to predict that a fraction of 
70% of all central star derives from a binary interaction. This frac- 
tion is even higher and more discrepant with the fraction detected 
here for the / band sample, but more in line with that derived from 
the J band sample. All these comparisons will need to be carried 
out again, once the PN binary fraction is finalised and a period dis- 
tribution is obtained. 



ACKNOWLEDGMENTS 

OD and J-CP acknowledge funding from NSF grant AST-06071 1 1 
which was used during the initial phases of this work. J-CP is 
thankful of Falk Herwig's support. We thank Thomas Rauch for his 
help in obtaining the best synthetic stellar atmospheres. We thank 
Dimitri Douchin for expediting his analysis of the variability prop- 
erties of the stars in our sample. We are grateful to Howard Bond 
for helping in the observing proposal writing stage. The TheoSSA 
service (http:dc.g-vo. orgtheossa) used to retrieve theoretical spec- 
tra for this paper was constructed as part of the activities of the 
German Astrophysical Virtual Observatory. 

REFERENCES 

Abell G. O., 1966, ApJ , 144, 259 



Afsar M., Bond H. E., 2005, Memorie della Societa Astronomica 
Italiana, 76, 608 

Bannister N. P., Barstow M. A., Holberg J. B., Bruhweiler F. C, 

2003, MNRAS ,341,477 
Benedict G. F, McArthur B. E., Napiwotzki R., Harrison T. E., 

Harris H. C, Nelan E., Bond H. E., Patterson R. J., Ciardullo R., 

2009, AJ, 138, 1969 
Bentley A. F, 1989, in Torres-Peimbert S., ed., Planetary Nebulae 

Vol. 131 of IAU Symposium, A search for cool companions of 

planetary nebula nuclei., p. 312 
Bergeron P., Wesemael F, Beauchamp A., 1995, PASP , 107, 1047 
Bessel M. S., 1990, Astrophysics & Space Science, 83, 357 
Bessell M. S., 1979, PASP , 91, 589 
Bessell M. S., 1991, AJ, 101, 662 
Bessell M. S., Brett J. M., 1988, PASP , 100, 1 134 
Bilfkova J., Chu Y.-H., Gruendl R. A., Su K. Y. L., De Marco O., 

2012, ApJ Supp. , 200, 3 
Blocker T, 1995, A&A , 299, 755 

Bond H. E., 1994, in Shafter A. W., ed., Interacting Binary Stars 
Vol. 56 of Astronomical Society of the Pacific Conference Se- 
ries, Binary Nuclei of Planetary Nebulae, p. 179 

Bond H. E., 2000, in ASP Conf. Ser. 199: Asymmetrical Planetary 
Nebulae II: From Origins to Microstructures Binarity of Central 
Stars of Planetary Nebulae, p. 1 15 

Bond H. E., Meakes M. G., Liebert J. W., Renzini A., 1993, in 
Weinberger R., Acker A., eds, Planetary Nebulae Vol. 155 of 
IAU Symposium, HST Observations of the Nuclei of EGB 6 
(0950+ 139) and Abell 58 (V605 Aql). p. 499 

Cardelli J. A., Clayton G. C, Mathis J. S., 1989, ApJ , 345, 245 

Ciardullo R., Bond H. E., Sipior M.S., Fullton L. K., Zhang C.-Y, 
Schaefer K. G., 1999, AJ, 118, 488 

Cudworth K. M„ 1973, PASP , 85, 401 

De Marco O., 2009, PASP , 121, 316 

De Marco O., Farihi J., Nordhaus J., 2009, Journal of Physics 

Conference Series, 172, 012031 
De Marco O., Hillwig T. C, Smith A. J., 2008, AJ, 136, 323 
De Marco O., Passy J.-C, Moe M., Herwig F, Mac Low M.-M., 

Paxton B., 201 1, MNRAS , 411, 2277 
De Marco O., Soker N., 201 1, PASP , 123, 402 
Dengel J., Hartl H., Weinberger R., 1980, A&A , 85, 356 
Dorman B., Rood R. T, O'Connell R. W„ 1993, ApJ , 419, 596 
Dreizler S., 1999, in R. E. Schielicke ed., Reviews in Modern 
Astronomy Vol. 12 of Reviews in Modern Astronomy, Spec- 
troscopy of Hot Hydrogen Deficient White Dwarfs, p. 255 
Duquennoy A., Mayor M., 1991, A&A , 248, 485 
Edgar R. G., Nordhaus J., Blackman E. G., Frank A., 2008, ApJL 
, 675, L101 

Frew D. J., 2008, PhD thesis, Department of Physics, Macquarie 
University. 

Frew D. J., Madsen G. J., O'Toole S. J., Parker Q. A., 2010, PASA, 
27, 203 

Frew D. J., Parker Q. A., 2007, in Asymmetrical Planetary Nebu- 
lae IV, IAC Elec. Pub. Do post-common envelope objects form a 
distinct subset of PNe? 

Frew D. J., Parker Q. A., 2010, PASA, 27, 129 

Frew D. J., Parker Q. A., Boijcic I., 2012, MNRAS , accepted 

Fulbright M. S., Liebert J., 1993, ApJ , 410, 275 

Garcfa-Segura G., Langer N., Rozyczka M., Franco J., 1999, ApJ 
,517,767 

Garcfa-Segura G., Lopez J. A., Franco J., 2005, ApJ , 618, 919 
Gesicki K., Zijlstra A. A., 2007, A&A , 467, L29 



18 O. De Marco et al. 



Gianninas A., Bergeron P., Dupuis J., Ruiz M. T., 2010, ApJ , 720, 
581 

Harris H. C, Dahn C. C, Canzian B., Guetter H. H., Leggett S. K., 
Levine S. E., Luginbuhl C. B., Monet A. K. B., Monet D. G, Pier 
J. R., Stone R. C, Tilleman T., Vrba F. J., Walker R. L., 2007, 
AJ, 133,631 

Hartl H., Dengel J., Weinberger R., 1983, Mitteilungen der As- 

tronomischen Gesellschaft Hamburg, 60, 325 
Henry T. J., McCarthy Jr. D. W., 1993, AJ, 106, 773 
Holberg J. B., 2009, Journal of Physics Conference Series, 172, 

012022 

Holberg J. B., Bergeron P., 2006, AJ, 132, 1221 
Holberg J. B., Magargal K., 2005, in Koester D., Moehler S., eds, 
14th European Workshop on White Dwarfs Vol. 334 of Astro- 
nomical Society of the Pacific Conference Series, Finding the 
Cool Companions of the PG DA White Dwarfs, p. 419 
Jacoby G. H, Ferland G. J., Korista K. T., 2001, ApJ , 560, 272 
Kaler J. B., Shaw R. A., Kwitter K. B., 1990, ApJ , 359, 392 
Kirkpatrick J. D., McCarthy Jr. D. W., 1994, AJ, 107, 333 
Koester D., Schonberner D., 1986, A&A , 154, 125 
Kraus A. L., Hillenbrand L. A., 2007, AJ, 134, 2340 
Kun M., 1998, ApJ Supp. , 1 15, 59 
Landolt A. U., 1992, AJ, 104, 340 
Landolt A. U., 2009, AJ, 137, 4186 

Lawrence A., Warren S. J., Almaini O., Edge A. C, Hambly et al. 

2007, MNRAS , 379, 1599 
Leggett S. K., 1992, ApJ Supp. , 82, 351 

Liebert J., Bergeron P., Holberg J. B., 2005, ApJ Supp. , 156, 47 

Liebert J., Green R., Bond H. E., Holberg J. B., Wesemael E, 
Fleming T. A., Kidder K, 1989, ApJ , 346, 251 

Liu X.-W., Barlow M. J., Blades J. C, Osmer S., Clegg R. E. S., 
1995, MNRAS , 276, 167 

Lopez J. A., Richer M. G, Garcfa-Dfaz M. T., Clark D. M., 
Meaburn J., Riesgo H, Steffen W., Lloyd M., 2012, Rev. Mexi- 
cana de Astronorma y Astrofisica, 48, 3 

Lupton R., Ivezic Z., 2005, in Seidelmann P. K., Monet A. K. B., 
eds, Astrometry in the Age of the Next Generation of Large Tele- 
scopes Vol. 338 of Astronomical Society of the Pacific Confer- 
ence Series, Experience with SDSS: the Promise and Perils of 
Large Surveys, p. 151 

Madappatt N., De Marco O., Nordhaus J., Wardle M., 2011, in 
Schuh S., Drechsel H, Heber U., eds, American Institute of 
Physics Conference Series Vol. 1331 of American Institute of 
Physics Conference Series, Towards a New Prescription for the 
Tidal Capture of Planets, Brown Dwarfs and Stellar Compan- 
ions, pp 319-323 

Mafz-Apellaniz J., 2004, PASP , 116, 859 

Mastrodemos N., Morris M., 1999, ApJ , 523, 357 

Miszalski B., Acker A., Moffat A. F. J., Parker Q. A., Udalski A., 
2009, A&A , 496, 813 

Miszalski B., Acker A., Parker Q. A., Boffin H. M. J., Frew D. J., 
Mikolajewska J., Moffat A. F. J., Napiwotzki R., 201 1, in Asym- 
metric Planetary Nebulae 5 Conference Dense circumstellar neb- 
ulae in wide binary central stars, p. 109P 

Moe M., De Marco O., 2006, ApJ , 650, 916 

Moe M., De Marco O., 2011, in Asymmetric Planetary Nebulae 
5 Conference Population synthesis of galactic PN from strong 
binary interactions, p. 107 

Morales-Rueda L., Maxted P. F. L., Marsh T. R., North R. C, 
Heber U., 2003, MNRAS , 338, 752 

Napiwotzki R., 1999, A&A , 350, 101 

Napiwotzki R., Schonberner D., 1995, A&A , 301, 545 



Nordhaus J., Blackman E. G„ Frank A., 2007, MNRAS , 376, 599 

Nordhaus J., Spiegel D. S., Ibgui L., Goodman J., Burrows A., 
2010, MNRAS , 408, 631 

Parker Q. A., Acker A., Frew D. J., Hartley M., Peyaud A. E. J., 
Ochsenbein E, Phillipps S., Russeil D., Beaulieu S. E, Cohen 
M., Koppen J., Miszalski B., Morgan D. H, Morris R. A. H, 
Pierce M. J., Vaughan A. E., 2006, MNRAS , 373, 79 

Passy J.-C, De Marco O., Fryer C. L., Herwig E, Diehl S., Oishi 
J. S., Mac Low M.-M., Bryan G. L., Rockefeller G., 2012, ApJ , 
744, 52 

Raghavan D., McAlister H. A., Henry T. J., Latham D. W., Marcy 
G. W., Mason B. D., Gies D. R., White R. J., ten Brummelaar 
T. A., 2010, ApJ Supp. , 190, 1 

Rauch T, Deetjen J. L., 2003, in I. Hubeny, D. Mihalas, & 
K. Werner ed., Stellar Atmosphere Modeling Vol. 288 of As- 
tronomical Society of the Pacific Conference Series, Handling 
of Atomic Data, pp 103 — h 

Rauch T, Koppen J., Napiwotzki R., Werner K, 1999, A&A , 
347, 169 

Rauch T, Werner K., 1995, in D. Koester & K. Werner ed., White 
Dwarfs Vol. 443 of Lecture Notes in Physics, Berlin Springer 
Verlag, New Analyses of Helium-rich Pre- White Dwarfs, p. 186 

Rauch T, Ziegler M., Werner K, Kruk J. W., Oliveira CM., 
Vande Putte D., Mignani R. P., Kerber E, 2007, A&A , 470, 317 

Rodriguez M., Corradi R. L. M., Mampaso A., 2001, A&A , 377, 
1042 

Rosado M., Moreno M. A., 1991, Astrophysics & Space Science, 
88, 245 

Schlafly E. E, Finkbeiner D. P., 2011, ApJ , 737, 103 
Schmidt-Kaler T, 1982, Bulletin dTnformation du Centre de 

Donnees Stellaires, 23, 2 
Schonberner D., 1983, ApJ , 272, 708 

Schonberner D., 1993, in R. Weinberger & A. Acker ed., Plane- 
tary Nebulae Vol. 155 of IAU Symposium, Evolutionary Tracks, 
p. 415 

Soker N, 1997, ApJ Supp. , 1 12, 487 
Soker N, 2006, PASP ,118, 260 

Tody D., 1986, in Crawford D. L., ed., Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series Vol. 627 of 
Society of Photo-Optical Instrumentation Engineers (SPIE) Con- 
ference Series, The IRAF Data Reduction and Analysis System, 
p. 733 

Tody D., 1993, in Hanisch R. J., Brissenden R. J. V, Barnes 
J., eds, Astronomical Data Analysis Software and Systems II 
Vol. 52 of Astronomical Society of the Pacific Conference Se- 
ries, IRAF in the Nineties, p. 173 
Tremblay P.-E., Bergeron P., Gianninas A., 2011, ApJ , 730, 128 
Tweedy R. W., Kwitter K. B., 1996, ApJ Supp. , 107, 255 
Vacca W. D., Garmany C. D., Shull J. M., 1996, ApJ , 460, 914 
Vassiliadis E., Wood P. R., 1994, ApJ Supp. , 92, 125 
Villaver E., Livio M., 2009, ApJL , 705, L81 
Wegner W., 1994, MNRAS , 270, 229 

Weinberger R., Dengel J., Hartl H, Sabbadin E, 1983, ApJ , 265, 
249 

Werner K., 1995, Baltic Astronomy, 4, 340 

Werner K., Deetjen J. L., Dreizler S., Nagel T, Rauch T, Schuh 
S. L., 2003, in I. Hubeny, D. Mihalas, & K. Werner ed., Stel- 
lar Atmosphere Modeling Vol. 288 of Astronomical Society of 
the Pacific Conference Series, Model Photospheres with Accel- 
erated Lambda Iteration, pp 31-+ 

Werner K., Dreizler S., 1999, Journal of Computational and Ap- 
plied Mathematics, 109, 65 



The Binary Fraction of Planetary Nebula Central Stars I. 



19 



Werner K., Herwig R, 2006, PASP , 118, 183 

Zuckerman B., Becklin E. E., McLean I, S., 1991, in Elston 
R., ed., Astronomical Society of the Pacific Conference Series 
Vol. 14 of Astronomical Society of the Pacific Conference Se- 
ries, Central stars of planetary nebulae in the infrared., pp 161- 
166 



APPENDIX A: INDIVIDUAL PHOTOMETRIC 
MAGNITUDES OF THE CENTRAL STAR OF PN SAMPLE 

In Table I A 1 1 we report all the photometric measurements for our 
central stars of PN. The measurements for each star were averaged 
according to the Equations[3]and[4]in §|4]and presented in Table|6] 



APPENDIX B: THE THEORETICAL JOHNSON-COUSINS 
COLOURS OF HOT STARS 

In this appendix we report the colours determined for a grid 
of stellar atmosphere models with Solar and PG1159 abun- 
dances. The theoretical stellar atmosphere models were calculated 
with the simulation code TMAW, the web i n terface to TMAP 
dWerner & Dreizler! 1 19991 : IWerneretal] 120031 : iRauch & DeetierJ 
l2003h . or the German Astrophysical Virtual Observatory grid cal- 
culations TheoSS^f] which are carried out for standard, non- 
tailored abundances but typically with a more complete model 
atom. TheoSSA is deemed more correct but differences were typ- 
ically lower than 1%. The TMAW stellar atmosphere models are 
calculated for a solar composition (with mass fractions /3_ff=0.71, 
/? ffe =0.28 /?c=0.001 /3iv=0.002 /3 o =0.006). The TheoSSA mod- 
els are for /?h=0.7, /3f/ e =0.3 unless indicated o r for the composi- 
tion t ypical of hydrogen-deficient PG1 159 stars dWerner & Herwi j 
120061 : /3 We =0.33 /3 C =0.50 /3 N =0.02 /3 o =0.15). The colours are 
then determined using the routine calcphot within the iraf pack- 
age (using the vegamag and effstim options). The filter system used 
to obtain the synthetic colours is the Johnson-Cousins system, the 
same used for the determination of the standard magnitudes by 
Landol 3 dl992h . These colours were modified by the zero points 
listed in table 3.1 of the synphot manual, 0. 010 for the B — V and 
R-I colours and -0.002 for the V - 1 one dMarz-Apellaniz_2004l ; 
iHolberg &~B ergeron 2006). The V — J and J — H colours have 
not been corrected. We note that the offset for the V band alone 
is 0.026 mag. If the correction for the J magnitude were zero, this 
would imply that the theoretical colours are redder, and the result- 
ing excess smaller by about 3%. We note that none of our detec- 
tions in the V — J band would be affected, although the detected 
spectral types would be slightly fainter. The theoretical colours for 
Solar as well as PG1 159 abundances are tabulated in Table IB li as a 
function of temperature and gravity, where we have also tabulated 
blackbody values for comparison. The parameter space covered is 
that which applies to central star of PN as well as post-RGB, hor- 
izontal branch sta rs. See for instance the log g — T c g diagram in 
lNapiwotzkiUl999l) . 



APPENDIX C: THE MAGNITUDES, COLOURS AND 
MASSES OF MAIN SEQUENCE STARS 

In this Appendix we report the value for absolute V magnitudes 
and colours for main sequence stars of spectral type B3V to LOV 
These data were used to interpret the I and J band excesses we 
detected in terms of companion spectral type. 

For the colours of the companions (assumed to be main se- 
quence stars) we have used a variety of data from the literature. For 
spectral types earlier than B8V, we av eraged the 1^-band ab solute 
magnitudes an d B — V colours from ISchmidt-Kalerl q 1983) and 
Wegner] dl994l) . with the data from Eric Mamajek's websitfQ (last 
updated in Dec 2011). For spectral types in the ra nge B8V-K7V, 
we av erage d the data from Mamajek's website, ISchmidt-Kaled 
i 19821) and iKraus & H illenbrand! j2007l) . transformed fol lowin 



zing 
ssel 



iLupton & Ivezid ( 20050, with additional B — V data from lBesse 
jl990h andlBessellldl99ll). We averaged the My and B — V data 
from lLeggetd jl992l)~ Kirkpatric k & McCarthy! dl994h . Mamajek's 
website and Frew (unpublished) for spectral typ es later than MO V. 

The V — R c and colours were taken from lBesse ll ( 1979) or 
iBessell d 199dh for spectral types earlier than MOV, and from lLeggettl 
dl992h . iKirkpatrick & McCarthy! d 19941) and Frew (unpublished) 
for spectral types l ater than MOV; The V — I c co l ours were simi- 
larly a v eraged fromlBe ssell ( 1979), Bessel ( 1990), Bes sell & Brettl 
dl988h . iLeggettl d 19921) , IKirkpatrick & McCarthy! dl994l) . Mama- 
jek's website, and Frew (u npublished). The V — J an d V — H 
colours were derived from Rra us & Hillenbrand! d2007t) and Ma- 
majek's website for spectral t ypes earlier than MOV, and from 
Kirkp atrick & McCarthy! ( fl994l) . Mamajek's website and this work 
for spectral types later or equal to MO V. Finally, the stellar masses 
were derived from the relations given in iHenrv & McCarthy! fl 99 3) 
for spectral types later than A5V, together with other literature es- 
timates for earlier spectral types. 



3 dc.zah.uni-heidelberg.de/theossa/ 



4 http://www.pas. rochester.edu/~emamajek/EEM_dwarf_UB VIJHK_colors_Teff.dat 

5 http://www.sdss.org/dr5/algorithms/sdssUBVRITransform.html 
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18.129 


18.173 


18.114 


18.052 


0.021 


0.023 


0.034 


0.055 


1.2 




4 


18.124 


18.168 


18.120 


18.101 


0.024 


0.015 


0.020 


0.025 


1.2 




5 


18.177 


18.200 


18.160 


18.092 


0.013 


0.012 


0.020 


0.028 


1.2 




6 


18.157 


18.156 


18.101 


18.139 


0.022 


0.023 


0.029 


0.050 


1.4 




7 


18.115 


18.113 


18.088 


18.102 


0.013 


0.011 


0.020 


0.025 


1.3 


JnEr 1 


4 


16.770 


17.143 


17.288 


17.483 


0.018 


0.014 


0.018 


0.020 


1.1 




7 


16.778 


17.117 


17.287 


17.531 


0.014 


0.010 


0.018 


0.033 


1.2 


K 1-13 


3 


18.062 


18.425 


18.569 


18.798 


0.013 


0.017 


0.021 


0.027 


1.3 




4 


18.045 


18.420 


18.603 


18.896 


0.013 


0.013 


0.020 


0.030 


1.2 




5 


18.036 


18.435 


18.613 


18.852 


0.033 


0.026 


0.037 


0.061 


1.2 


K2-2 


4 


13.984 


14.278 


14.396 


14.536 


0.011 


0.013 


0.017 


0.020 


1.1 




5 


13.974 


14.261 


14.395 


14.561 


0.012 


0.010 


0.020 


0.025 


1.4 




7 


13.972 


14.254 


14.379 


14.564 


0.011 


0.010 


0.018 


0.019 


1.3 


NGC 3587 


4 


15.414 


15.771 


15.953 


16.164 


0.013 


0.014 


0.019 


0.020 


1.3 




5 


15.415 


15.788 


15.963 


16.227 


0.014 


0.012 


0.018 


0.027 


1.3 




6 


15.413 


15.768 


15.966 


16.207 


0.021 


0.023 


0.027 


0.047 


1.3 


NGC 6720 


4 


15.393 


15.793 


15.898 


16.050 


0.012 


0.020 


0.020 


0.019 


1.2 




7 


15.420 


15.747 


15.904 


16.075 


0.016 


0.019 


0.022 


0.022 


1.6 


NGC 6853 


4 


13.725 


14.075 


14.244 


14.404 


0.012 


0.016 


0.019 


0.020 


1.3 




5 


13.777 


14.100 


14.241 


14.390 


0.012 


0.013 


0.019 


0.045 


1.1 




7 


13.744 


14.090 


14.255 


14.417 


0.013 


0.013 


0.020 


0.035 


1.6 


PuWe 1 


5 


15.298 


15.551 


15.673 


15.802 


0.013 


0.010 


0.019 


0.023 


1.4 




7 


15.285 


15.539 


15.651 


15.785 


0.011 


0.010 


0.018 


0.018 


1.3 


Sh 2-78 


4 


17.643 


17.666 


17.630 


17.567 


0.014 


0.016 


0.030 


0.038 


1.4 




5 


17.624 


17.655 


17.580 


17.513 


0.013 


0.012 


0.038 


0.046 


1.3 


Sh 2-176 


3 


18.639 


18.580 


18.597 


18.602 


0.022 


0.024 


0.036 


0.059 


1.1 




4 


18.460 


18.569 


18.571 


18.521 


0.018 


0.017 


0.020 


0.029 


1.1 




7 


18.431 


18.536 


18.553 


18.541 


0.011 


0.015 


0.020 


0.028 


1.1 


Sh 2-188 


4 


17.438 


17.454 


17.407 


17.361 


0.012 


0.014 


0.017 


0.019 


1.1 




5 


17.416 


17.446 


17.402 


17.381 


0.012 


0.009 


0.018 


0.025 


1.2 




7 


17.419 


17.444 


17.385 


17.386 


0.013 


0.009 


0.017 


0.017 


1.1 


Ton 320 


5 


15.385 


15.731 


15.900 


16.126 


0.013 


0.011 


0.020 


0.025 


1.3 




7 


15.373 


15.718 


15.880 


16.091 


0.015 


0.012 


0.019 


0.037 


1.2 


WeDe 1 


5 


16.961 


17.229 


17.329 


17.503 


0.018 


0.011 


0.023 


0.028 


1.7 




6 


16.962 


17.221 


17.334 


17.470 


0.012 


0.011 


0.018 


0.027 


1.3 




7 


16.946 


17.231 


17.354 


17.500 


0.021 


0.022 


0.027 


0.046 


1.1 



Table Al. The individual photometric magnitudes for our sample of PN central stars (including three PN mimics). The date and time corresponding to each 
night can be found in Tabled DeHt 5 and HDW 4, having only one epoch of data are reported only in Table|6] 
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T eff 


log g 


B — V 


V — I 


V — J 


R — I 


J — H 


Abundance 


(kK) 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




20 


4 


-0.195 


-0.193 


-0.466 


-0.097 


-0.073 


Solar (TMAW) 


20 


5 


-0.180 


-0.199 


-0.475 


-0.099 


-0.076 


Solar (TMAW) 


20 




-0.139 


-0.118 


-0.414 


-0.066 


-0.068 


BB 


30 


4 


-0.269 


-0.277 


-0.653 


-0.145 


-0.109 


Solar (TMAW) 


30 


5 


-0.260 


-0.280 


-0.660 


-0.145 


-0. 1 1 1 


Solar (TMAW) 


30 


7 


-0.221 


-0.289 


-0.680 


-0.143 


-0.106 


Solar (TMAW) 


30 


8 


-0.191 


-0.289 


-0.695 


-0.139 


-0. 1 00 


Solar (TMAW) 


30 


~ 


-0.226 


-0.214 


-0.582 


-0.116 


-0.098 


BB 


40 


4 


-0.271 


-0.281 


-0.684 


-0.150 


-0.114 


H0.728He0.249 


40 


5 


-0.279 


-0.298 


-0.713 


-0.158 


-0.119 


H0.894He0.106 


40 


6 


-0.274 


-0.304 


-0.724 


-0.161 


-0.121 


H0.738He0.014 


40 


7 


-0.277 


-0.319 


-0.759 


-0.165 


-0.126 


Solar (TMAW) 


40 


8 


-0.263 


-0.318 


-0.759 


-0.161 


-0.121 


Solar (TMAW) 


40 




-0.266 


-0.259 


-0.662 


-0.140 


-0.113 


BB 


50 


4 


-0.280 


-0.291 


-0.704 


-0.156 


-0.116 


H0.738He0.249 


50 


5 


-0.290 


-0.305 


-0.725 


-0.162 


-0.119 


H1.0 


50 


6 


-0.298 


-0.310 


-0.745 


-0.164 


-0.125 


Solar (TMAW) 


50 


7 


-0.294 


-0.316 


-0.758 


-0.166 


-0.127 


Solar (TMAW) 


50 


8 


-0.280 


-0.315 


-0.760 


-0.166 


-0.126 


H0.7He0.3 


50 


~~ 


-0.289 


-0.285 


-0.708 


-0.154 


-0.122 


BB 


60 


5 


-0.302 


-0.315 


0.753 


-0.169 


-0.126 


H0.7He0.3 


60 


6 


-0.307 


-0.319 


-0.755 


-0.169 


-0.129 


Solar (TMAW) 


60 


7 


-0.305 


-0.322 


-0.770 


-0.170 


-0.130 


Solar (TMAW) 


60 


8 


-0.295 


-0.322 


-0.771 


-0.170 


-0.129 


H0.7He0.3 


60 




-0.304 


-0.302 


-0.737 


-0.163 


-0.128 


BB 


70 


5 


-0.310 


-0.325 


-0.771 


-0.174 


-0.130 


H0.7He0.3 


70 


6 


-0.312 


-0.323 


-0.722 


-0.171 


-0.131 


Solar (TMAW) 


70 


7 


-0.311 


-0.328 


-0.780 


-0.173 


-0.132 


Solar (TMAW) 


70 


8 


-0.306 


-0.326 


-0.783 


-0.175 


-0.132 


H0.6He0.4 


70 




-0.314 


-0.314 


-0.760 


-0.169 


-0.132 


BB 


80 


5 


-0.316 


-0.331 


-0.781 


-0.177 


-0. 1 30 


H0.7He0.3 


80 


6 


-0.318 


-0.332 


-0.787 


-0.178 


-0.133 


H0.7He0.3 


80 


7 


-0.318 


-0.333 


-0.791 


-0.176 


-0.135 


Solar (TMAW) 


80 


8 


-0.314 


-0.334 


-0.797 


-0.179 


-0.135 


H0.7He0.3 


80 




-0.322 


-0.323 


-0.774 


-0.174 


-0.135 


BB 


90 


5 


-0.316 


-0.330 


-0.779 


-0.177 


-0.129 


H0.7He0.3 


90 


6 


-0.326 


-0.340 


-0.806 


-0.181 


-0.140 


Solar (TMAW) 


90 


7 


-0.325 


-0.340 


-0.804 


-0.180 


-0.138 


Solar (TMAW) 


90 


8 


-0.322 


-0.341 


-0.808 


-0.183 


-0.137 


H0.7He0.3 


90 




-0.328 


-0.329 


-0.786 


-0.177 


-0.137 


BB 


100 


6 


-0.332 


-0.347 


-0.817 


-0.184 


-0.141 


Solar (TMAW) 


100 


7 


-0.331 


-0.346 


-0.814 


-0.183 


-0.139 


Solar (TMAW) 


100 


8 


-0.329 


-0.345 


-0.814 


-0. 1 85 


-0.138 


H0.7He0.3 


100 




-0.332 


-0.335 


-0.796 


-0.180 


-0.139 


BB 


110 


6 


-0.330 


-0.341 


-0.801 


-0. 1 82 


-0.132 


H0.7He0.3 


110 


7 


-0.336 


-0.349 


-0.819 


-0.185 


-0. 1 40 


Solar (TMAW) 


110 


8 


-0.334 


-0.347 


-0.815 


-0. 1 86 


-0.138 


H0.7He0.3 


110 


~ 


-0.336 


-0.339 


-0.803 


-0.182 


-0.140 


BB 


120 


6 


-0.332 


-0.343 


-0.803 


-0.183 


-0.132 


H0.7He0.3 


120 


7 


-0.339 


-0.352 


-0.822 


-0.186 


-0.141 


Solar (TMAW) 


120 


7.5 


-0.337 


-0.347 


-0.812 


-0.186 


-0.136 


H0.7He0.3 


120 


~~ 


-0.339 


-0.343 


-0.810 


-0.184 


-0.142 


BB 


130 


6 


-0.333 


-0.344 


-0.804 


-0.183 


-0.132 


H0.7He0.3 


130 


7 


-0.341 


-0.353 


-0.827 


-0.187 


-0.142 


Solar (TMAW) 


130 


8 


-0.340 


-0.350 


-0.818 


-0. 1 87 


-0.137 


H0.7He0.3 


130 




-0.342 


-0.346 


-0.815 


-0.186 


-0.143 


BB 


140 


6 


-0.335 


-0.345 


-0.805 


-0.184 


-0.132 


H0.7He0.3 


140 


7 


-0.342 


-0.354 


-0.830 


-0. 1 88 


-0. 143 


CiOJai ( 1 MAW ) 


140 


8 


-0.343 


-0.351 


-0.819 


-0.188 


-0.137 


H0.6He0.4 


140 




-0.344 


-0.348 


-0.820 


-0. 1 87 


-0. 144 


BB 


150 


6 


-0.335 


-0.346 


-0.808 


-0.184 


-0.132 


HO. /He0.3 


150 


7 


-0.344 


-0.355 


-0.831 


-0.188 


-0.143 


Solar (TMAW) 


150 


8 


-0.343 


-0.352 


-0.820 


-0.188 


-0.137 


H0.8He0.2 


160 


6 


-0.336 


-0.347 


-0.8 10 


-0. 185 


-0. 132 


HO. /HeO. 5 


160 


7 


-0.345 


-0.356 


-0.833 


-0. 1 89 


-0. 144 


CiOJai ( 1 MAW ) 


160 




-0.345 


-0.353 


-0.820 


-0.188 


-0.136 


H0.7He0.3 


160 


6 


-0.335 


-0.349 


-0.814 


-0.185 


-0.134 


H0.8He0.2 


170 


7 


-0.346 


-0.358 


-0.837 


-0.190 


-0.145 


Solar (TMAW) 


170 


8 


-0.346 


-0.353 


-0.821 


-0.189 


-0.136 


H0.7He0.3 


170 




-0.349 


-0.354 


-0.830 


-0.190 


-0.146 


BB 


100 


6 


-0.307 


-0.321 


-0.769 


-0.178 


-0.151 


He0.33C0.50N0.02O0.15 


100 


7 


-0.340 


-0.357 


-0.835 


-0.194 


-0.151 


He0.33C0.50N0.02O0.15 


110 


6 


-0.369 


-0.35 1 1 


-0.987 


-0.231 


-0.204 


He0.33C0.50N0.02O0.15 


110 


7 


-0.346 


-0.355 


-0.832 


-0.194 


-0.148 


He0.33C0.50N0.02O0.15 


120 


6 


-0.386 


-0.382 


-0.876 


-0.203 


-0.143 


He0.33C0.50N0.02O0.15 


120 


7 


-0.350 


-0.357 


-0.836 


-0.195 


-0.148 


He0.33C0.50N0.02O0.15 


130 


6 


-0.387 


-0.384 


-0.878 


-0.204 


-0.145 


He0.33C0.50N0.02O0.15 


130 


7 


-0.352 


-0.359 


-0.838 


-0.196 


-0.147 


He0.33C0.50N0.02O0.15 



This colour was interpolated between the values for the 1 lOkK and 130kK atmospheres. 



Table Bl. Predicted colours of single post-AGB stars using TMAP models (TMAW with solar metallicity - see text, or TheoSSA, where the model atom and 
abundances are indicated as mass fractions adjacent to the relevant atom) as well as for blackbody curves (BB). These values have been adjusted for the small 
corrections of 0.010 mag for the B — V and R — I colours and -0.002 mag for the V — I colours 
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Spec. M v (U-B) (B-V)o (V - R c )o (V - I c )o (V - J)o (V - H) (V - K 3 ) Q Mass 
type (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) M 



B2 V 


-2.10 


-0.82 


-0.23 


-0.09 


-0.23 


-0.49 


-0.57 


-0.60 


7.7 


B3 V 


-1.45 


-0.69 


-0.19 


-0.07 


-0.19 


-0.38 


-0.47 


-0.49 


5.7 


B4 V 


-1.22 


-0.63 


-0.18 


-0.06 


-0.18 


-0.34 


-0.43 


-0.45 


5.1 


B5 V 


-1.06 


-0.58 


-0.16 


-0.05 


-0.17 


-0.32 


-0.41 


-0.42 


4.8 


B6 V 


-0.77 


-0.50 


-0.15 


-0.04 


-0.15 


-0.28 


-0.35 


-0.36 


4.2 


B7 V 


-0.53 


-0.44 


-0.13 


-0.04 


-0.11 


-0.25 


-0.32 


-0.33 


4.0 


B8 V 


-0.24 


-0.35 


-0.10 


-0.03 


-0.09 


-0.19 


-0.26 


-0.26 


3.5 


B9 V 


+0.16 


-0.20 


-0.07 


-0.02 


-0.06 


-0.09 


-0.14 


-0.12 


3.1 


AO V 


+0.79 


-0.02 


-0.01 


+0.00 


-0.01 


+0.05 


+0.01 


+0.04 


2.4 


Al V 


+1.05 


+0.03 


+0.03 


+0.01 


+0.03 


+0.09 


+0.07 


+0.10 


2.2 


A2 V 


+1.36 


+0.06 


+0.06 


+0.03 


+0.07 


+0.17 


+0.16 


+0.19 


2.1 


A3 V 


+1.53 


+0.08 


+0.09 


+0.05 


+0.10 


+0.20 


+0.20 


+0.23 


2.0 


A4 V 


+1.74 


+0.09 


+0.13 


+0.06 


+0.17 


+0.30 


+0.32 


+0.36 


1.95 


A5 V 


+1.90 


+0.10 


+0.15 


+0.07 


+0.19 


+0.36 


+0.37 


+0.40 


1.89 


A7 V 


+2.16 


+0.10 


+0.20 


+0.12 


+0.24 


+0.45 


+0.49 


+0.53 


1.76 


A8 V 


+2.35 


+0.09 


+0.25 


+0.15 


+0.29 


+0.51 


+0.59 


+0.63 


1.65 


FOV 


+2.63 


+0.04 


+0.30 


+0.18 


+0.36 


+0.60 


+0.69 


+0.73 


1.58 


F2 V 


+3.00 


-0.00 


+0.36 


+0.21 


+0.43 


+0.72 


+0.88 


+0.93 


1.46 


F5 V 


+3.46 


-0.02 


+0.43 


+0.27 


+0.52 


+0.89 


+1.03 


+1.08 


1.31 


F8 V 


+4.01 


+0.01 


+0.53 


+0.30 


+0.60 


+1.06 


+1.23 


+1.29 


1.19 


GOV 


+4.40 


+0.05 


+0.58 


+0.33 


+0.67 


+1.12 


+1.36 


+1.42 


1.10 


G2 V 


+4.72 


+0.12 


+0.64 


+0.36 


+0.71 


+1.18 


+1.47 


+1.55 


1.03 


G5 V 


+5.07 


+0.19 


+0.67 


+0.39 


+0.74 


+1.25 


+1.56 


+1.64 


0.99 


G8 V 


+5.51 


+0.29 


+0.74 


+0.43 


+0.78 


+1.38 


+1.69 


+1.77 


0.93 


KO V 


+5.89 


+0.44 


+0.82 


+0.46 


+0.85 


+1.50 


+1.86 


+2.03 


0.87 


Kl V 


+6.08 


+0.52 


+0.85 


+0.50 


+0.88 


+1.55 


+1.94 


+1.95 


0.85 


K2 V 


+6.37 


+0.62 


+0.90 


+0.54 


+0.93 


+1.67 


+2.06 


+2.16 


0.82 


K3 V 


+6.61 


+0.80 


+0.98 


+0.57 


+1.05 


+1.81 


+2.30 


+2.41 


0.80 


K5 V 


+7.34 


+1.08 


+1.16 


+0.73 


+1.29 


+2.19 


+2.71 


+2.84 


0.72 


K7 V 


+8.16 


+1.20 


+1.34 


+0.81 


+1.57 


+2.54 


+3.24 


+3.41 


0.65 


MOV 


+8.87 


+1.19 


+1.41 


+0.89 


+1.76 


+2.86 


+3.50 


+3.74 


0.59 


Ml V 


+9.56 


+1.18 


+1.47 


+0.96 


+1.98 


+3.20 


+3.77 


+4.03 


0.54 


M2 V 


+10.17 


+1.17 


+1.50 


+1.00 


+2.14 


+3.36 


+3.94 


+4.18 


0.45 


M3 V 


+11.01 


+1.17 


+1.55 


+1.08 


+2.45 


+3.80 


+4.38 


+4.62 


0.33 


M4 V 


+12.80 


+1.18 


+1.67 


+1.19 


+2.75 


+4.41 


+4.96 


+5.23 


0.24 


M5 V 


+14.20 


+1.3 : 


+1.82 


+1.41 


+3.30 


+5.13 


+5.73 


+6.00 


0.15 


M6 V 


+16.59 


+1.3 : 


+2.03 


+1.81 


+3.93 


+6.25 


+6.86 


+7.19 


0.11 


M7 V 


+17.84 




+2.15 


+2.13 


+4.51 


+7.03 


+7.64 


+8.02 


0.10 


M8 V 


+18.72 




+2.15 


+2.24 


+4.57 


+7.55 


+8.23 


+8.67 


0.09 


M9 V 


+19.39 




+2.15 


+2.37 


+4.61 


+7.72 


+8.45 


+8.94 


0.08 


LOV 


+19.65 






+2.44 


+4.66 


+8.10 


+8.82 


+9.36 


0.08 



Table CI. Spectral types, empirical V-band absolute magnitudes, intrinsic colour indices and masses of main-sequence stars 



